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SYNOPSIS 


In the present investigation the influence of powder 
injection parameters, for example mass flow rate of the powdered 
flux, carrier gas flow rate etc. » on the kinetics of impuricy 
extraction via a mathematical model is studied. 

Effect of impurities of interest i . e. > sulphur, phosphorous 
and nitrogen > on steel properties is presented in chapter 1. .4 

mathematical model is developed for two cases; that is with slag 
carry-over and without slag carry-over A general solution is 
obtained for the case 'with slag carry-over". Whereas FDM is used 
for the case "without slag carry-over" for the reasons given in 
chapter 2. I nvesti gati on on model parameters selection is done. 
These are also presented in Chapter 2. 

Chapter 3 is reserved for the study of model parameters on 
concentr ati on vs time curves under the two secuions 'with slag 
carry over ^ and "without slag carry-over", transitoric race 
contribution in the overall refining process and its dependency/ cn 
these model parameters and transitoric efficiency. . 

In chapter 4 applications of developed model are presentee. 
It involves the evaluation of model parameters from the given 
concentr at i on vs time curves ^ and from the exper i mentals resulis 
collected from the literature by using the selected equations in 
chapter 3. 

Conclusions and suggestions are given in chapier 5 and 
chapter 6 respectively/. Solutions procedures of developed model 
equations with the slag carry/-over and without slag carry-over are 
presented in Appendix —1 and Appendix -2 respectively. Appendix - 
3 gives detailed description of transitoric rate contribution. 
Programme code for model results is presented in Appendix - 5 and 
the procedure is documented in flow-charts 1 and 2. 
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CHAPTER 1 
INTRODUCTION 

A tremendous increase in the interest and development of 
refining processes by injection metallurgy, including powder and 
gas injection, has found growing acceptance since it calls for 
quicker processing at reasonable reaction efficiency for a wider 
span of heat size. This is to meet the demands for high quality 
products like pipe line steels. Applications of injection 
metallurgy are mostly found in steel processing units. Therefore, 
a brief review on physical metallurgy of impurities, with respect 
to sulphur, phosphorus and nitrogen in steel is required, 

1.1 Physical Metallurgy of Impurities CS^ N, 1 in Steels: 

Increasing demand of steel materials for more stringent 
applications have called for much efforts to decrease impurity 
elements to extra-low concentrations. For example extra-deep 
drawing high tensile strength strip, super-ferritic stainless 
steel strip, hydrogen-induced-crack (HIC) and 
sulphide-stress— corros ion— crack (SSCC) resistant high tensile 
strength plate and pipe, Ci — Mo pressure vessel plate, 9^. Ni plate 
for cryogenic temperature service etc. 

Generally impurities present either as dissolved in the 
matrix, often not uniformly distributed but enrich to certain 
structural features like grain boundaries, or in the form of NMI - 
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1.1.1 Sulph.urs 

Except in free machining steels, sulphur is the highly 
undesirable element in steels. It has been so well known that 
sulphur causes hot-shortness when the steel is forged or rolled at 
elevated temperature by forming low melting point eutectic alloy 
i.e., FeS at the grain boundaries. It also deteriorates the 
ductility of steel materials via precipitation during 
solidification of and deformation during rolling of manganese 
sulphides and pitting corrosion in Austenitic stainless steels. 

Among the properties relevant to represent the ductility, 
shelf energy at low temperature and reduction of area <RA) in 
through thickness direction are most markedly influenced by the 
sulphur. For example in heavy reactor pressure vessel plates used 
in nuclear power plant, the shelf energy in V-notch charpy test 
improves from 19 to 22 Kgf-m with the decrease of sulphur from fcO 
to 30 ppm Cll. The values of RA in through- th ickness direction of 
lamellar tear resistant carbon steel plate linearly go up as 
sulphur reduces from 50 to 20 ppm when dissolved oxygen falls in 
between 10 to 20 ppm C21. 

The most important problems with sulphur in hydrogen 
environment is the formation of HIC and SSCC. Formation of HIC 
and SSCC have been well known in HSLA steels used for pipe lines.. 
Effective measures to prevent their formation is to control both 
the sulphur concentration in the melt and the shape of sulphide 


inclusions 
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SSCC is also caused by the pitting corrosion often observed 
in Austenitic stainless steels. Origin of pitting corrosion is 
located to be NMI among which manganese sulphides are most 
harmful. It was confirmed that pitting potential of 18*/. Cr - 8*/. 
Ni steel increases when steel composition is changed from 1.0*/. hn 
and 50 ppm S to 0,2V, Mn and i 10 ppm S C33. 

Thus for high quality steel sulphur should be kept as low as 
possible . 

1.1.2 Phosphorus* 

Next to sulphur, phosphorus is the undesirable element in 
steels. It has the tendency to segregate to grain boundaries. It 
has been so well known to deteriorate ductility of steel material, 
thereby increasing the tendency of steel when cold worked, making 
it cold-short. It also known to cause temper embrittlement, HIC, 
Stress-Corrosion-Cracking (SCO etc. 


Steel materials most 

sensitive 

to 

phosphorus 

is 9% 

Ni steel 

plate used for 

cryogenic 

temperature 

serv ice . 

Even 

very 

low 

concentrations 

of phosphorus make 

9*/. 

Ni steel 

susceptible 

to 


temper embrittlement and deteriorate low temperature ductility. 
Even for air cooled specimen more than 5 Kgf-m absorbed energy at 
~196°C is gained by decreasing phosphorus from 130 to 80 ppm. 
Similarly crack opening displacement <C0D) also improve from 0.06 
to 0.2 as phosphorus decreases from 80 to 10 ppm, for the furnace 
cooled specimen C23. Temper embrittlement of low carbon Ni-Cr-Mo 
steels will be appreciably deteriorated by phosphorus. Forged 
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plates of these steels show improvement in ductile to brittle 
transition temperature is improved from 100®C to 0°C at 0.3*/* Si 
and 50 to -75*^0 at 0.02X Si by decreasing phosphorus from 150 to 
20 ppm C43. Similarly by decreasing phosphorus from 170 to 60 ppm 
notch toughness of reactor pressure vessel plate in nuclear power 
plant can be improved from -150 to -175°C. In all these grades of 
steels decrease of phosphorus decreases intergranular segregation 
and hence improves their granular embrittlement, contributing to 
increased ductility of steels. 

Similar to sulphur phosphorus also causes HIC of HSLA steels 
used for line pipes to transport sour gas and air. In case of 
sulphur the culprit was shape of sulphides, whereas in case of 
phosphorus local enrichment of phosphorus. If the steels are to 
be free from HIC, phosphorus must be less than 0.1*/* for tin of 
2.2*/. The same threshold values are applicable to prevent the 
formation of SCC in this grade of steel C23. 

Transgranul ar SCC susceptibility of the austenitic stainless 
steels is extremely improved by decreasing phosphorus from lOO to 
20 ppm. 

Besides low temperature properties of steels, it also 
deteriorates hot ductility of steels especially when the 
concentration is above lOO ppm. The embrittlement at P > 100 ppm 
is due to the formation of liquid film or brittle eutectic 
phosphide at austenitic grain boundaries. The remedy is to 
decrease phosphorus below 100 ppm and/or cooling the strand cast 
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slabs in the secondary cooling zone at slower rate to avoid 
segregation of phosphorus. 

Therefore decrease of phosphorus to very low levels (tenth 
order ppm) is essential to produce high quality steels. 

1.1.3 Nitrogens 

For steel products like deep drawn steels and line pipes very 
low concentration of nitrogen in steel is required, since affects 
yield point elongation, ductility, and corrosion of steels. 

Decrease in nitrogen is effective to prevent yield point 
phenomenon, a drop of elongation to fracture. When nitrogen is 
decreased from 80 to 20 ppm the sharp increase in the ductility of 
weld bonds of Al- killed steel plate was observed E21. 

Presence of nitrogen (and carbon) causes intergranular 
corrosion of super ferritic stainless steel strips. Their 
limiting value decrease as the concentration of Cr increases. For 
example 130 ppm of Ni+C for 26*/. Cr and 80 ppm for 30*/. Cr steel. 
Such mal-inf luence of nitrogen <and carbon) is due to the 
formation of Cr-depleted zone around the intergranularly 
precipitated chromium carboni trades. 

The decrease of nitrogen to extra-low concentrations helps 
considerably improving yield point elongation, elongation to 
fracture and pitting corrosion of steels. 

The demand for the production of high quality steel materials 
changed the philosophy of steelmaking. Powder injection is being 
increasingly employed by the steel industry since it involves less 
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investment, is easy to operate and facilitates quicker processing 
at reasonable reaction efficiency for a wide span of heat sizes. 
Despite its growing acceptance, the kinetics of refining processes 
are poorly understood. The kinetics have a direct impact on the 
quality of the product and the cost of production, because the 
particles reside in the melt for very less time say of the order 
of I sec or so. Thus the identification of process variables is 
necessary to study the kinetics of powder injection refining. 

1.2 Literature Survey 

Usually sulphur and phosphorus are removed from steels via 
slag/metal reactions, whereas Nitrogen removal is done by 
degassing. But there have been several reports of Nitrogen 
removal from steel by the addition of flux and to regard nitrogen 
removal as one of the slag/metal reactions. Slag/metal reactions 
of interest are desulphurisation (hereafter termed as DeS) , 

dephosphorisation (hereafter, termed as DeP> and deni troni zat ion 
(hereafter termed as DeN> • 

The efficiency of slag/metal reactions depends on how good 

the metal and slag are in contact with each other. The following 

are the types of contact modes of slag with metal. 

(i) Permanent contact 

(ii) Transitoric contact 

(iii) Counter current contact 

Schenck C53 defined the removal without relative movement of 
the phases as permanent phase contact and with relative movement 
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of the phases as transitoric phase contact or counter current flow 
process* In transitoric phase contact only one phase is mobile, 
whereas in counter current flow process both phases move in 
opposite direction to each other. 

These three contact modes are schematically shown in Fig, 
1.1, Fig. 1.2 and Fig. 1.3. The contact modes of interest are <i> 
and < i i ) on 1 y . 

As the production of high quality steels require decrease of 
impurity to low concentrations, powder injection refining is being 
widely employed in steel industry. In the plant as well as 
laboratory scale sufficient experimental work was done in the 
powder injection, refining. In the literature powder injection 
refining was purely correlated to the transitoric contact mode. 
But practically the top slag formed from the injected flux is not 
removed. Hence its influence on the powder injection refining is 
to be considered. Only Robertson et al . £63 take into account the 
top slag contribution in the powder injection refining process via 
permanent contact reaction. Although these are two more reference 
(£73, £83) on powder injection taking into account the top slag 
influence on the kinetics of refining, they are based on Robertson 
et al , work . 

In spite of the development of powder injection refining 
process, it is still unexplained whether permanent contact mode or 
transitoric contact mode dominates the process. Although 
Robertson et al., tried to describe the powder injection refining 
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via a mathematical model, the effect of injection variables were 
not described explicitly. The solution method they specified is 
FDM. This model requires the prefixing of permanent rate constant 
i.e., Kp value, Robertson et al . , assumed a Kp value to fit their 
experimental data in the model. In the Ref. C73 and C83 Kp was 
prefixed by the permanent contact mode experiments. 

In Robertson et al . , model the experimental points can be 
fitted by many number of combinations of Kp and E. <Kp and E are 
the two variable in their model). Thus when Kp is assumed to be 
certain value, some value of E will be obtained by fitting the 
experimental data points in the model. There are also some works 
where Kp value is fixed by assuming weak agitation and strong 
agitation C93. Despite the uncertainty in the determination of 
critical agitation ( or stirring energy) assumption of weak 
agitation and strong agitation and there by fixing the Kp value 
involve more uncertainty in the validity of the results. 

In the process of metal holding in a ladle resulphurisation 
and rephosphorisat ion, i.e, reversion of sulphur and phosphorus 
from slag to steel are reported C93,C103. 

The sources of resuphurisation and rephosphorisat ion are 

(a) carry-over slag from the previous refining unit 

(b) erosion of ladle lining 

The effect of carry-over slag on slag/metal refining is more 
prominant then the erosion of ladle lining because of the 
improvement in the quality of ladle lining refractories. Although 
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use of devices such as refractory sphere, slide gate etc, and EBT 
(Eccentric Bottom tapping) in EAF control slag carry-over in the 
ladle, but can not eliminate it completely. E.T. Turkdogan 
discussed the reactions that occur between liquid steel and slag 
in the ladle C103. But the discussion is confined to theoretical 
considerations only. He presented the effect of furnace slag 
carry-over on the recovery of A1 added and rephosphorisation in 
the tap ladle. In the ref C93 resulphurisat ion due to FeO contact 
of the overall slag is presented. 

The experimental conditions of powder injection as well as 
dumping (or permanent contact mode) are compiled from literature 
and are presented in Table 1.1 and Table 1.2 for hot-metal and 
steel respectively. The data from their tables is used in the 
following chapter to obtain the range of variation of powder 
injection parameters like loading ratio, stirring energy etc. 

1.3 Aim of the Present work* 

The present work is aimed at 

(i) The development of a mathematical model to obtain melt 
concentration as a function of time. 

(ii) The study of various parameters on transitoric mode of 
contact contribution and the conditions under which permanent 
or transitoric contact mode dominate the refining process. 
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CHAPTER 2 

MODEL DESCRIPTION 

In this chapter the basis for model formulation and the 
numerical solution procedure for the governing differential 
equations with slag carry-over from BQF and without slag 
carry-over, the significance of model parameters and the selection 
of equations are presented. 

2.1 Basis for Model Formulation: 

In slag/metal refining processes required amount of slag 
forming materials are either added to the top surface of the 
molten metal or injected directly into the bath at some injection 
rate. Out of these two types of additions, powder injection is 
being increasingly employed by the steel industry since it 
facilitates quicker procesing of wide span of heat sizes and high 
recovery of vaporizing fluxes Like CaC^, CaSi etc. Usual slag 
forming materials are mixtures of either CaO, CaSi^? ^^ 2 ^ 3 ’ ^^^2 
or calcium based compounds such as CaC^, CaSi , CaAl or magnesium 
and or sodium carbonate and lime or calcium carbonate and lime. 
Inspite of its increasing applications, the kinetics of powder 
injection refining processes are poorly understood. The kinetics 
have a direct impact on the quality of the product and the cost of 
production because the injected particles reside in the melt for 
very less time say 1 sec or so. Thus to maximize the extent of 
refining, refining rate and the reagent utilization, the effect of 
operating variables such as solid flow rate, gas flow rate, type 
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of fluK etc., must ba clearly understood. 

Thus the present model is developed to study the kinetics of 
slag/metal reactions caused by the injection of slag forming 
materials with carrier gas via a lance submerged into the molten 
bath. The gas-particle mixture penetrates into the melt until 
their momentum is dissipated. The powder particles detach from 
the gas and disperse into the bath, and the particles start to 
rise to the top surface of the bath because they are buoyant. The 
carrier gas due to byoyancy effect accelerates rapidly. The drag 
force between gas and liquid causes the stirring in the bath. The 
accumulated injected flux on the top surface of the melt is not 
removed in powder injection processes. All these events are shown 
in Fig. 2,1. The physico-chemical characteristics of the 
resultant top slag would depend on whether there is slag present 
on the top of molten metal prior to start of the powder injection 
or not. Such slag may be carried over from the previous 
steelmaking unit eg. BOF . In the absence of slag carry-over the 
characteristics of the resulting top slag and injected flux are 
same. The accumulated top slag also takes part in the refining 
process. 

The removal of an impurity from metal can occur during the 
rise of slag particles and by the slag which is accumulating on 
the bath surface. The former mode of impurity removal is 
transitoric whereas the later is permanent mode. Thus the total 
rate of impurity removal in case of powder injection can be 



Carrier Gas + Powdered Flux 





Schematic model of" powje^ injection refinC^i 
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expressed by 

R = + Rp <2. 1 ) 

Where R^ is rate of removal of impurity during the rise of 
powder particles and Rp is the rate due to accumulated slag» 

Usually the slag-metal reaction rate is expressed fay the mass 
transport controlled first order rate equation 



Equation <2.2> contain two mass transport steps in series, 
one from the metal to slag interface, and the other into the bulk 
of the slag from the interface. 


Although the ratio of (3 to ( 3 - is 5-10, use of high partition 

w 

coefficient slags (rjg = 100-20001 in refining processes makes 

* 

much smaller than 1/^^. Thus we can make the assumption 
that metal side mass transfer resistance is the rate controlling 


step. Therefore the equation (2.2) can be rewritten by using 


dCxl 

dt 


ft ^ icxl 


■4o J 


(2.3) 


This is the general form of mass transfer equation applicable 
to the rate of removal of impurity whenever slag is contacted with 
the metal. Using equations (2.1) and (2.3) the following mass 
transfer equation in case of powder injection can be derived. 
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dC 

d 


-I . k, - . 4^ L-3 - ^’1 

ml J «« I J 


(2.4> 


where Cx ' 3 » Cx3 


< x^ > 


(2.5> 


is the interfacial area for transitoric reaction and can 
be determined by powder size and powder injection rate because the 
components of the powder are assumed to be non-volatile and do not 
dissolve in the melt during its entire refining period. 

Ap is the interfacial area for the permanent reaction which 
depends on bath stirring intensity produced by the carrier gas 
when it leave the melt. When emulsification of top slag into the 
melt occurs, Ap depends on the mass of slag emulsified and the 
drop size. In the absence of emulsification the interfacial area 
available for the reaction is simply the cross-sectional area of 
the melt container or the ladle. 

Concentration of impurity in the top slag at any instant of 
time is given by 

< 2 . 6 ) 

Combining equations <2. 4), (2.5) and <2. 6), we get 


(x) = (x^) + (Xp) 


dCx3 

dt 


'P - i'l . - i|'} 

^ I y)c; J "'I J 


(2.7) 


The mass balance equation are 

dCx3 


mg(xT) = 


m 


L dt 


( 2 . 8 ) 
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Wgt ( X ) 




tx ] - Cx3 
o 


(2.9) 


Equation (2.8) is the mass balance equation for the rising 
flux. Equation (2.9) is the mass balance equation for the 
resulting top slag and implies that the decrease in total impurity 
concentration in the melt is equal to increase in impurity 
concentrat ion in the accumulated slag. 

By putting equations (2.0) and (2.9) into (2.7), we get 


dCx3 

dt 


^ f 

I 


txl + 


dCxl 


* m„ dt 

— — 

”l 


1 “ 


( 2 . 10 ) 


Defining the variables 




K, 


®P = 


/5pAp 

__ , 

m 

^p^p 

Vm 


K, 


t. B_ = 
in j ’ T 




t. . 
in j 


T = 


t Cx3 

: 1 « * 


in 3 


Cx 3 
o 


and substituting them into equation (2.10), we get 


m, t . 

L in a 


( 2 . 11 ) 


dot 

dr 


= ot 


Pt Bp Cl -.XgT-J Bp/Xg 

1 ®T 1 f ®T ' 

t 1 +-^‘f J ll+ i f 


where 




( 2 . 12 ) 


(2. 13) 
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is the transitoric rate contribution in 
injection refining process. 


Defining 


®t ^ 


B, 



and ©g 





the 


overal i 


and substituting them into equation <2.12) » we get 


powder 


<2.14) 


1. - 1, I * !|| = - <2. IS) 

Equation <2.15) is the governing differential equation 
describing the decrease of impurity with time in powder injection 
refining without slag carry-over from EOF. 

3.1.2 With Carry-over Slag: 

In most cases there occurs some carry-over of slag while 
transferring molten steel from the steelmaking unit into the ladle 
in which further powder injection refining is planned to be 
performed. The carry-over slag is normally oxidizing in nature 
and is a mixture of CaO, SiO^t Al^Og, FeO, MgO etc. 
When powdered flux is injected into the molten baths on its rise 
to the top surface of the melt it will be mixed into the already 
existing top slag. As both top slag and powdered flux are of 
different composition the overall top slag composition changes 
through-out the injection period. Thus the partition coefficient 
of the resultant top slag is different from that of the injected 
flux. Indicating top slag partition coefficient as equation 
<2. 7) can be modified as 



IB 


dCx3 

dt 


^ f - - ^>1 


(2.16) 


The masws balance equations are 

• , , dCx3 

mgCXj) = - 


(rngt+m^,) (x > 

Defining Y_ = — 
u (n. 


r 1 

= m, Ax = m, Cx 3 - Cx3l 
L L o I 


( 2 » 8 > 

(2. 17) 

(2. 18) 


Putting equation (2»8) , (2ell)> <2.17) and (2.18) in (2.16), we 

get 


dot 

dr 


= ot 


;B +B fl - -;j 5 

1 1 < 


B, 




Substituting <2.14) into (2,19), we get 


d <<x- 


dT 




T 




1 + 


!=! 

vj 


I 




<«-l ) = 


<2,19) 


(2.20) 


Equation (2.20) is the governing differential equation for 
concentration change of the impurity in the melt with time during 
the powder injection when slag is carried-over from BOF prior to 
start of the injection. In the absence of slag carry-over Y^.= 0 
and T 7 ^ = thus the equation ’ <2.20 ) wi 11 be reduced to equation 
<2.15). 

To solve differential equations (2.15) and (2.20) the 


following assumption is must i.e.. 


the melt concentration is 
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uniform in the bulk metal. 

2. 2 Solution Procerdures: 

' In equation (2,15) time dependent variable is only a, where 

as in equation (2.20) there are two time dependent variables i.e, 

ot and 7}^. Analytical solution for equation (2.15) can be 

obtained. But the presence of two time dependent variables in 

equation (2,20) forced to adapt numerical solution method. 

2.2.1 Without Slag Carry-overs 

The form of differential equation (2.15) is very much similar 

to that of Leibnitz's linear equation 

^ + Yf (x) = C (2.21 ) 

dx 

Standard solution for equation (2.21) is available and can be 
obtained from any standard mathematics book. Due to their 
similarity in the form of differential equation^ solution of 
equation (2.21) can be used to get the general solution for 
equation (2.15). The detailed solution is given in Appendix 1» 
and the final solution is 


a . 1 - (^gT) exp (- .p^T) ) „,(1 ^ e ■> „,) 

m=0 


( 2 . 22 ) 


The summation series in equation (2.22) is convergent type 
and the maximum number of terms required would depend on the 
magnitude of as per the following equation 


N In 


[-^ 1 - 


In (N + 1)! = - 7.601 (2.23) 

CN = maximum number of terms! 


upto the third decimal accuracy. 



20 


2 , the maximum number of terms 

to be considered in the summation series of equation (2.21) are 9. 
2 . S. 2 With Slag Carry-overs 

To integrate the equation (2.20), equation of 17 ^ as a 
function of time must be available. But it is impossible to get 
single correlation of as a function of time for all types of 
slags. But it is possible to evaluate in time steps with the 
knowledge of mass flow rate of flux, composition of injected flux 
and top flux composition at the beginning of injection process. 
The suitable solution procedure is FDM. The available methods to 
solve equation ( 2 . 20 ) are 

(a) Implicit FDM 

(b) Explicit FDM 

(c ) Crank-Nicholson 's FDM 


For example when 


I 1 

(.1 + —Li 


Out of these three solution procedures, Crank-Nicholson ' s FDM 
is chosen to solve equation (2.20). Detailed solution is given in 
Appendix 2. 

The final solution is 




21 


Where 


T = 


At 


4'c 


t <>» > 3 _ , 

o av j+1 


(T ) . , + 

av j+l Y- 


(T > 

av 3+1 


T + T . 

3 + 1 3 


(y)*) 

av 


ifi*) ...+ (y>*) . 

o 3+1 3 


3 + 1 


; At = Time Step 


<2. 25) 


2. 3 Approxiroalions Front General Solution: 

Approximations, under the given conditions, are necessary to 
avoid unnecessary calculations. Once general solution is obtained 
it is easy to make approximations from it. 


2. 3.1 


B, 


00 


— ‘inj 

m 


CO 


As V and t. . are finite, the extreme condition B- — ► oo, 
is equivalent i.e. injection of fine powders at higher mass flow 
rate. By putting this condition in equation (2.14) one can get <ft- 


= Xg and Gg = O and the equation <2.22) will be reduced to 
a = exp (- XgT> 

2. 3. 2 X g > 00 


(2.26) 


i.e,, either — ► oo or Y — ► oo 

Thus the condition can be understood in terms of use of 
either higher partition coefficient flux or large amount of flux. 
By adding this condition to equation (2.14), 4^^ * + Bp, ©g= 0. 
Hence equation (2.22) will be approximated to 
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ot = exp 


(Bp + T 


(2.27> 


a. 3. 3 Eqiiillbrium Concentration at each ts 

(i) At equilibrium, Rate = 0 i.e,, = O 

aT 

(ii) Equilibrium can be achieved through the resulting top slag 
only i.e., at each t equilibrium concentration can be 
calculated fay putting = 0 

From equation (2.12) by substituting equilibrium condition 
1 


a. 


eq 


1+ \gT 


(T < 1) 


If T > I, T » 1 
t 


iX 


eq 


1 + X- 


( 2 . 28 ) 


(2.29) 


a. 4 Selection and Significance of Parameters: 


From equation (2.24) 

« 


a = ot (<)!) , e. 


S' 




(2.30) 


From equation (2.14) and (2.11) 



Aj, 


Ap, 




t . . , 
in 3 ’ 

f ) 

(2.31) 

= 8g(/?^, 

1 

/?p. 

Ap, 


* 


t . . , 
in 3 ’ 

f ) 

(2.32) 


T = /?(t . . ) 

in j 


(2.33) 

In Appendix - 3 detailed description of f is given. 

Dependency of transitoric rate contribution in the overall 


process of refining is given by 
f * f(/3^. A^, ^p, Ap, V^, Yg) 


(2.31) 


Usually t. Y , Y_, V are specific to a system. Thus the 
inj’ c’ S’ m ' 

parameters to be evaluated are r^g, /5^. A^, /?p, Ap. In the 

following sections selection of equation for each parameter as a 
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function of process variables is done. 

2.4.1 Partition Coefficient of the Slagy'Metal Systems 

Partition coefficient characterizes the state of the dynamic 
equilibrium with respect to the transferring entity. At 

equilibrium partition coefficient attains a maximum value and 
would depend upon the composition of slag and melt, and the 
operating temperature of the system. The transferring entities 
under consideration in the present study are sulphur, phosphorus 
and nitrogen. Hence the partition coefficient for the above 
entities are given below. 

Partition coefficient of the impurity is defined as the ratio 
of solute <or transferring entity) concentration in the slag to 
that in the melt. 

T) - j- ^ 2 (2.35) 

2. 4. 1 . 1 Sulphur: 

Sulphur distribution between molten metal and slag can be 
represented by an electrochemical reaction 

S + <0^~) <S^“) + g <2. 36) 


Equation (2.36) can be written as the combination of 

I <s^y + <0^~> {S^“> + I 

I is^y = §. 

0 = I <02> 

For equation (2.37) sulphide capacity is defined as 




(2.37) 

(2.38) 

(2.39) 


(wty. S) 


(2.40) 
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‘^<2.34) " '•^82 

*^<2.35) "^^02 

Putting equations (2.41) and C2.42) in (2.40) 


(2.41 ) 

(2.42) 


(2.43) 

(2.44) 

(2.45) 

(2.46) 
equation 

(2.47) 
Basic 

slags have high Cg. Thus from equation (2.47) the favourable 
conditions for deS are drawn to be 
— > high slag basicity 
— > low a^ or reducing atmosphere 

Due to low concentration of solute elements or dissolved 
elements in the mild steel grade fg can be taken as 1 , where as in 
hot metal it is evaluated to be nearly 3. 

In the literature several correlations are given for sulphide 
capacity evaluation as a function of slag composition and 
temperature. All of them are listed in Table 2.1. 

Equation (2) in Table 2,1 employs optical basicity concept. 


Cg= (wty,S) , K ^ 2 . 34 >‘^(2,35) 

From the definition of partition coefficient > 7 , 

Combining equations (2.43) and (2.44) 
a 

* 2 - 

*"S “ %u 7^ '^(2.34)'^(2.35) 


Su 


(%S) 

cysi 


log = log Cg - log 


f J ’^(2.34)'^ '^(2.35)j 

j t_._ ?.^6 j I 


log K ^2 34 )'*' ^(2 35) given by - 1.375 and thus 

(2.46) becomes 




log ■S 7 g^ = log Cg - log 


936 


+ 1.375 


C^ is strong function of composition of slag system. 
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Optical basicity of slag system can be evaluated by the fallowing 
equation: 



where X * Mole fraction of oxide or fluoride 
afbt..c>d....f Oxidation numbers 

= Optical basicity of oxide or fluoride given in Table 

2 , 2 . 


Activity of oxygen is controlled by the type and amount of 

deoxidizer. To reduce soluble oxygen to very low level and to 

have low deoxidizer conccontration in the melt a strong deoxidizer 

must be added. Normally A1 is used to reduce the level of oxygen. 

Thus in equation <2.47) a can be correlated to the A1 in the melt 

o 

by the reaction 


2 A1 + 30 e* <Al202> <2,49) 

log 45 ) = - 20.63 <2.50) 

In equation <2.49) it is assumed that the melt is in 

equilibrium with the suspended Al^O^ particles, i.e. a^^ q =1. By 

2 3 

putting the a^ evaluated from equation <2.50) in the equation 


<2.47) the following equation 


2 

log = log Cg + ^ log 




<2.51 ) 


can be obtained 
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6x10 


In case of carbon saturated iron a^ is approximately equal to 
-5 


The optical basicity correlation (Equation (2) in Table 2.1> 
is used in this work to evaluate sulphide capacity and thus 
sulphur partition coefficient by equation (2.51). 
a. 4. 1 . 2 Phosphoruss 

Phosphorus distribution between molten metal and slaq can be 
represented by an electrochemical reaction 

P + I Q + I (0^“> = <P0^") (2.52) 


Phosphorus capacity for the above reaction (2.52) is usually 
given by 


<y.P) 

P " 5/2 


(2.53) 


From the definition of partition coefficient 

•/?p = (‘/,P)/C‘/«P1 (2.54) 

Combining equation (2.53) and (2,54) 

* 5 

loQ yip - log *-p 2 ^ 12.55 ) 

log fp = 0 for mild steels 

= 0.6043 for hot metal 

Same as sulphide capacity, phosphorus capacity is also strong 
function of slag system and is high for highly basic slags. Thus 
from equation (2.55) the favourable conditions for deP are drawn 
to be 


: high slag basicity 

! high a or oxidizing atmosphere 

Similar to sulphide capacity, phosphorus capacity 
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correlations are also available in the literature, as a function 
of slag composition and temperature , and are listed in table 2.3. 
In correlation <2> of Table 2.3 optical basicity can be evaluated 
as explained in section (2.4. 1.1). Correlation 1 of Table 2.3 
gives activity coefficient of PO^ ^ formed from the reaction 


P + 2.50 = (P0„ _) 

— ~ 2.5 

= 11 ^- 8.51 

Phosphorus capacity can be obtained by the 


procedure s 


K 




(2.56) 2.5 

a . a 

B. a 


2.5 

a„ . a 
P o 


(•/.P ) 


PQ 


2.5 


2.5 MoHei 

P Slag 


(2.56) 

(2.57) 
f ol lowing 


It^ Cp 



® log 


PO 


2.5 


' Vsiag s 
*^( 3 . 56 ) ” 



(2.58) 


where Mp = atomic wt, of Phosphorus 

n , = number of moles of slag 

slag ^ 

Optical basicity correlation (Equation (2), Table 2.3) is 
used in this work to evaluate Cp , But for slag systems containing 
special components, for example CaCl^, the regular solution model 
with the interaction energies of the components collected from 
literature is used in Cp calculation. 
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2. 4. 1 . 3 Nitrogen: 

Nitrogen distributes between slag and raetal according to the 


electrochemical reaction 

N + I (O^") = (N^T + I Q. 

Equation <2.59) can be written as the combination of 


i + I <o^"> 

i tN2> * N 


(N^ ) + J 


AB^ = 3590 + 23.89T <J/mole) 

I = Q 

.o = -117150 - 2.887T )J/mole) 

2 . q 2 

For equation (2.60) nitride capacity is defined by 
p3/4 

3- °2 

V- = ^ ^ 

N„ 


(2.59) 

(2.60) 
(2,61 ) 

(2.62) 


(2,63) 


K 


(2,61 ) 


N 

t?F 


N, 


(2.64) 


K 


(2,62) 


o 

Wr 


(2,65) 


Putting (2.64) and (2.65) in (2.63) we get 

r -.3/2 


Cj^3- = (%N^ ) 


K 


£2 

a, 


r 0 1 ' 

.61 ) - 

N [’^(2.62) 


From the definition of partition coefficient 




(7,N) 


'N cy.N3 

Combining (2.67) and (2.66) , we get 


log = log C^Z- - ^ log a^+log f,^ + 5 log K^2.62)" ‘^(2.61) 




( 2 . 66 ) 


(2.67) 
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3 H 

= log C|^3- - ^ log + log ' ■ + 1.4742 (2.68) 

According to equation (2.66), it is evident that in order to 
obtain a high value of nitrogen distribution ratio» it is required 
that the slags have high nitride capacity and the metal low oxygen 
activity. In contrast to sulphide capacity, nitride capacity 
decreases with increasing basicity of slag. Nitride capacity for 
two slags systems is shown (Fig. 2.2 and 2.3), 

As explained in the section (2, 4.1.1) a^ can be correlated to 
the • Thus the fallowing equation relating partition 

coefficient with nitride capacity and dissolved aluminium 
activity 

log n* = log c„3- - I log log 4 log 8.8406 

(2.69) 

can be obtained. 

For some slag systems nitrogen partition coefficient as a 
function of dissolved Al content in the melt is shown in Fig. 2.4. 

In the nitrogen partition coefficient calculation, because of 
unavailability of correlations, nitride capacity values from Fig. 
2.2 and 2.3 are used. 

3. 4. S Tr-ansitoric rate constant CJCj.I>s 

To differentiate the rate constant terms appearing in 
equations 2.22 and 2.24, the one designated by K.^. i.e., rate 

constant for the rising flux, is termed as transitoric rate 
constant and the other designed by Kp i.e, rate constant for the 
accumulated top slag, is termed as permanent rate constant. 
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Fig. E.3 Nitride Capacity of CaO-BaO-Al 0 -TiO at 1873K 

S3 e 







Nitrogen partition 
sol . A1 



33 


Transitoric rate constant depends on transitoric mass transfer 
conefficient ifij) (similar to K^, ftj means the mass transfer 
coefficient for the rising flux), transitoric area (A^) i.e., 
slag/melt interfacial area of the rising flux and the volume of 
the melt. As the volume of the melt is specific to a system, the 
main parameters are fij and A^. Rate of impurity extraction by the 
rising flux depends on the magnitude of this rate constant besides 
the driving force term i-e., concentration gradient term in 
equation (2.4). Hence, the parameters on which depends, are 
discussed in this section. 

2. 4. 2.1 Transitoric mass transfer coefficient C^^>i 

Mass transfer coefficient is the furnction of diffusivity of 
the solute, fluid properties, injected powder-properties etc. Due 
to its complex dependency on these variables it is usually 
expressed by empirical correlations. The magnitude of transitoric 
mass transfer coefficient is strongly influence by the state 
of injected flux i.e., whether it is in solid or liquid form 
during its rise in the molten metal. Depending upon the melting 
point of the flux it might remain in the solid or liquid form 
during its rise. When melting point of the flux is more than the 
refining temperature i.e, 1623K in case of hot-metal and i873K in 
steel, injection of premelted flux is useful, since it melts at 
faster rate, upon its contact with the molten melt. 

When injected flux melts in the molten bath, it forms liquid 
droplets. There are large number of correlation available in the 
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literature to describe mass transfer in a wide variety of 
liquid-liquid system involving drops. These correlations are 
listed in Table 2.4 and compared in Fig. 2.5, A band is observed 
in this Figure with the equation 1 <table 2.4) for solid spheres 
as lower boundary and equation 5 (Table 2,4) for circulating and 
oscillating drops combined with penetration theory as upper 
boundary. Drops rising through the melt have high degree of 
internal circulation caused by stirring, and also as per 
penetration theory, there will be surface renewal at liquid-liquid 
interface caused by bubbles. That is why equation 5 from Table 
2.4 is used in this work for the rising drops in the molten metal. 

When injected flux is in the solid form on its contact with 
the melt, the solid reaction product formed at the flux/metal 
interface slows down the overall mass transfer coefficient. For 
the mass transfer between solid and liquid phase Asai et al, C543 
proposed that 


1 2 
^Sh « '^Re ^Sc 

where a^ and a^ are arbitrary constants 


(2.70) 


When this type of situation is encountered the equation for 
solid spheres i.e,, equation I from Table (2.3) is the probable 
solution in the absence of other data. 

Apart from fluid properties and diffusivity of the solute in 
the melt, the terminal velocity of the droplets in the melt must 
be known to evaluate It can be evaluated as per the following 


procedure . 
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2. 4. a. 2 Terndnal velocity of the droplets In the melt: 

When the drops reach the state of terminal velocity, the sum 
of all the forces acting on the drops i,e., the summation of 
gravity force, buoyancy force and force due to fluid motion, must 
be zero. 

Then , 


3 ^m^ 


4 _3 

3 ^s« 


u 


m 


(2,71 > 





To solve equation (2,72) for terminal velocity, C^j must be 
known. The approximate equation given for is C533. 

By solving equation (2.72) and (2,73) simultaneously for the 
given drop diameter and fluid properties, terminal velocity can be 
obtained . 


2. 4.2. 3 Trans! toric interfacial area CA^>: 

Transitoric interfacial area is the effective interfacial 
area available in the residence time t^ when the powdered flux of 
diameter dp is injected at the mass flow rate iHg. The magnitude 
of interfacial area is affected by the amount of particles wetted 
by the molten metal. In the calculation of A^, from the following 
correlation, it is assumed that the whole injected flux is 
available for the reaction. Therefore, 
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L 1 Kl 

^ - [dFj 


<2.74) 


Residence time of injected flux in the bath is the ratio of 
lance submergence depth to the velocity of the rising flux. 

Usually lance submergence depth, if not given, is 507. of the melt 

depth. Thus the only unknown is the velocity of the rising flux 

in the molten metal. The particles rise with terminal velocity 
with respect to liquid. As the carrier gas stirr the bath, the 

rising velocity of the particles can be obtained by adding its 

terminal velocity and the molten metal velocity. As the molten 

metal velocity is much greater than the terminal velocity ofthe 

particles, particle velocity can be taken as the molten metal 

velocity. To calculate the molten metal velocity, the equation 
(2.75) C553 derived for the mean rising velocity of the molten 
metal in the plume zone is used. 

Mean rising velocity of the molten metal in the plume zone 

^ 2 ^ I ■.O.b2 


Uj. = 19.9 


PI 


'R 


Ql 


0 


<2.75> 


<2.76) 


2.4.3 Permanent Rate Constant CKp>t 

Similar to transitoric rate constant permanent rate constant 
also depends on mass transfer of solute from bulk metal to the top 
slag, total interfacial area <i.e., ladle cross-sectional area and 
emulsified area) and the volume of the melt. As the equilibrium 
conditions put a constraint on the driving force term 
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effectiveness of top slag in the refining processes will be 
decided by the magnitude of the rate constant. That is higher the 
magnitude of the rate constant, faster will be the refining 
reactions. In a slag/metal system the main contribution to the 
permanent rate constant will be from mass transfer coefficient and 
total interfacial area. The entrainment of top slag into the 
molten metal enhances both and Ap . But unfortunately 
quantitative data on entrainment volume, number, size of droplets 
etc., is not available. In the literature empirical correlations, 
on mass transfer coefficient, for systems related to the 
(undisturbed) area of liquid/liquid interface and for systems 
involving drops (see Table 2.4) are available. But due to the 
uncertainity in the determination of drop diameter and entrained 
mass of top slag, rate constant is usually corelated to the 
stirring energy supplied via carrier gas to the molten metal. For 
this purpose the rate curves presented in the literature for 
various slag systems (see table 2.5) for permanent contact mode 
(alternatively termed as dumping) experiments were evaluated for 
Kp. The evaluated data show that Kp vary in between 0.0093 to 

0. 55 depending upon the type of slag system, state of the system 

1. e, whether it is in liquid or solid state, and the stirring 
intensity. Thus the evaluation of Kp as a function of stirring 
energy, from this table alone, will be erroneous. Therefore Kp is 
evaluated, in gas stirring, from the following two correlations. 

(a) for Phosphorus C563 
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Kp = 7,0x10“^ (2.77) 

<b) for Sulphur C571 

Kp = 7.6xl0~^<js;)^*^ (2.78) 

Equation (2.77) was derived by using oxidizing basic slags 
with £ values in the range lO to 500 W/T in 250T ladle, whereas 

m 

equation (2.78) was derived at e values in the range 60 to 100 W/T 
with 2.5T liquid steel. Depending upon the requirement 
extrapolation or interpolation of equations (2.77) and (2.78) is 
done . 

In inductively stirred melts, where no breakage of top 
slag/metal interfacial area occurs, /3p varies from 
0.002 - 0.003 Cm/sec. C381. 

Thus the correlations selected in this section for the 
calculation of model parameters are used in later chapter for the 
applications of model to study deS, deP and deN. 
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CHAPTER 3 

RESULTS AND DISCUSSION 

In chapter 2 the equations for powder injection refining 
process with and without slag carry-over were derived (Equation 
2.22 without slag carry-over and Equation 2.24 with slag 
carry-over). With slag carry-over and without slag carry-over are 
discussed explicitly to give their distinguishing features. This 
chapter is devoted to the parameters influencing change in melt 
concentration with time and also their effect on rate of impurity 
extraction by the top slag and the rising flux. 

3, 1 No Slag Carry-Over: 

From equation 2.22 it can be seen that the melt concentration 
with time depends on variables ^ and d. Furhter and depend 

on Bp» Xg and f. As and dg depend on the same parameters y 

they can not be varied independently. Thus Bp, Xg and f are 

the variables to be studied. As given in Appendix-3, f is 
function of B_, B_ and X_; thus the parameters influencing melt 

FI o 

concentration with time are reduced to Bp, B^, Xg. Xg can be 
varied fay varying composition of flux and/or mass fraction of flux 
to melt. Bp can be varied by varying time of injection and/or 
stirring energy i.e., carrier gas flow rate. B^ can be changed by 
changing time of injection and/or mass flow rate of flux and/or 
lance immersion depth and/or powder particle size. 

In refining processes metal sampling is done in time 

intervals to know whether the end point is reached or not. Hence 
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the concentration vs time curve arises from this process. To give 
salient features of concentration vs time curves, it is necessary 
to plot them as a function of independent variables on which they 
depend. And also it is necessary to know rate of change of melt 
concentration with time. Thus the effect of X_, B_and B_on 

o r I 

concentration change with time are discussed independently. 

Evaluation of concentration as a function of time for the given 

B_, B_ and X_ is shown in Flow-chart 1. 

P T S 

3.1.1 Influence of Flux Type and/or its amount: 

The influence of type of flux is characterized quantitatively 
by the partition coefficient. Both and Yg are grouped in Xg. 

The effect of Xg on dimensionless concentration change with time 
is shown in lower part of the Fig. <3. la). The upper part shows 
the variation of rate as a function of time for two selected 
values of Xg. The plot is divided into two periods, t from O to 

1 is the injection period and 1 onwards it is the stirring period. 

The selected Bp and B^ values to study Xg effect on a vs t curve 
are 0.4 and 0.9 respectively. Xg is varied between 1 to 18. 
Analysis of Xg effect above 18 is unnecessary because of 

negligible difference between concentration curves between Xg=15 
and Xg=18. This is the indication of ineffectiveness of Xgabove 
15. The equilibrium curves, evaluated from equation <2.28) and 
equation (2.29), are also superimposed in this plot. 

The increasing difference between equilibrium and kinetic 

curves with increasing Xg value can be observed from Fig. 3.1(a). 
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This indicates that Xg affects the contribution from the rising 
flux in the injection period^ otherwise the difference between 
equilibrium and kinetic curve would have been the same. And also 
the increase in difference in stirring period indicates that the 
top slag must be effective at higher Xg value in the overall 

refining process. This observation will be confirmed in the later 
section where Bp effect on concentration curve is discussed. The 
rapid decrease in concentration when Xg is increased from 1 to 5 
is due to increase in driving force for impurity removal. Above 
15 the increase in driving force is negligibly small. This is due 
to the negligible contribution from <x)/y> for the driving force 
term in the Eq. 2.4. The independent rate variation i.e., 

permanent rate, transitoric rate and total rate with time, is 
shown in Fig. 3.1<b> for two Xg values i.e., and Xg*5m At 

T=1 the difference between transitoric and permanent rates at 
is less than that at Xg=l. That is, the top slag is effective at 
higher Xg values. The increase in magnitude of rate indicates 
that increase in stirring might be helpful. This will be analysed 
in the following section. 

3.1.2 Infltiezice of Transitoric Parameter CB^>s 

B.J. involves both ^.nd KLj.. Prolonged time of refining is 

undesirable from kinetic point of view- Analysis of literataure 

experimental data suggests that nearly 10 min in powder 

injection refining processes. Hence the variation of B^ must be 
understood as the variation of Although involves volume of 
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the melt, mass transfer coefficient, mass-flow rate of the flux 
and residence time, its significant variation will be affected by- 
mass flow rate of the flux and residence time only. Hence its 
variation can be thought of varying mass flow rate of the flux 
and/or residence time of powder particles in the molten bath. The 
effect of transitoric parameter on dimensionless concentration as 
a function of dimensionless time is shown in Fig, 3.2<a> and Fig. 
3.3<a) for two X. values 1 and 5, 

From Fig. 3.2(a) and Fig. 3.3(a) three types of values can 
be identified, depending upon the melt concentration in relation 
to the equilibrium concentration at t=1. 

i. giving the melt concentration above the equilibrium 

concentration at t=1. 

(e.g.s B^=0,6, 0.9 in Fig. 3.2(a) and B^“0.6, 0.9, 1.9 

in Fig. 3, 3(a) ) . 

ii. B^ giving the melt concentration equal to the equilibrium 
concentration at t=1 

(e.g.: B^ = 0.19 in Fig. 3.2(a)>. 

iii. B^ giving the melt concentration below the equilibrium 

concentration at t=1 

(e.g.t B^=4.16 in Fig. 3.2(a> and Fig. 3.3(a)) 

Corresponding rate curves are given in Fig. 3.2(b) and Fig. 
3.3(b), 

For B^=4.16 in Fig. 3.2(b) and Fig. 3.2(b) negative rate can 
be observed i.e., impurity is transferring from slag to the melt. 




and its rate of removal as a function cf time 
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This is the phenomena of back transfer. Thus the first type 
i.e, B^=0,6, 0,9 in Fig, 3,2<a) and B^=0.6, 0,9, 1.9 in Fig, 
3.3<a> can be termed as ' B^ corresponds to no back-transfer 
condition', second B^ type i,e. B^®1.9 in Fig. 3.2<a> is ' B^ 
corresponds to equilibrium conditions' and termed as optimum B^, 
and the third B^ type i.e., B^=4.l<b in Fig. 3.2(a> and Fig. 3. 3(a) 
can be termed as ' B^ corresponds to back transfer condition'. In 
case of 'B^ corresponds to back transfer condition' any further 
stirring, after the injection period, is unnecessary because in 
this region top slag deteriorates the advantage of injection 
processes by way of back transfer. For 'B^ corresponds to 
equilibrium condition' also no further strirring is needed, 
because the top slag attains equilibrium value at t=: 1 and can not 
extract the impurity. Stirring is advantageous only, for ' B^ 
corresponds to no back transfer condition'. 

In Fig. 3.3(a), the departure of kinetic curve from 
equilibrium curve is observed to be larger than that in Fig. 
3.2<a) and the delay in the intersection of kinetic and 
equilibrium curve can be observed. Although there is large 
decrease in concentration value when Xg is increased from 1 to 5, 
but the increased difference between kinetic and equilibrium curve 
suggests that there might be decrease in transitoric effeciency. 
This also indicates that the top slag contribution can be 
increased at higher X„ values by way of stirring. 

Hence in injection process the effect of top slag must be 
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3.4 Effect of B on impurity concetration in the melt 

P 

and its rate of removal as a function of time (B_ =0.9) 
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than at Xg-1. Thus the large shift in concentration curve at X_®5 

S 

when Bp is varied from 0.4 to 1.8 at B^=0,9 indicates that 
stirring effect is more dominant at higher X^ values. A glance at 

w 

Fig. 3.€ gives the view that the top slag is not effective for 
B^ values giving excess impurity removal, in injection period, and 
after the injection period back transfer is proportional to the Bp 
value. Observation of Fig, 3,4<a) and Fig. 3.S causes to make 
the conclusion that Bp is effective parameter to reach equilibrium 
melt concentration. That is why top slag is only effective when 
the kinetic curve is above the equilibrium curve. This fact is 
also shown in Fig, 3.4(b> and Fib. 3.5 , where rate is plotted 

as a function of dimensionless time. In Fig. 3.4<b> i.e. 

B^=0.9 rising flux is still dominant as its rate at all values of 
T is above the rate of impurity extraction correspond to the top 
slag, whereas in Fig. 3.5.^, for Bp*1.0 top slag rate is higher 
than that of rising flux upto t is 0.4. Hence the higher stirring 
intensity, when mass flow rate of flux is low, makes the top slag 
as effective as the rising flux. 

3.1.4 Transitoric rate Contribution and Trans! toric Efficiency: 

In the literature there is no agreement on the domination of 
permanent contact reaction of transitory contact reaction in the 
powder injection refining process. Therefore rate contribution 
from either of them and the factors influencing the rate 
contribution, relation between transitoric efficiency parameters 
and the transitoric rate contribution have been included in this 
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chapter . 

3. 1.4. 1 Transltoric rate contributtom 

The method to ofatian transitoric rate contribution is given 

in flow-chart 1 for the given Bp, and values. It was 

already discussed that it depends on B_, B_ and only, B value 

r T S T 

is selected from no back transfer region, because above this B^ 

value, Bp has no effect in the injection period as given in 

section 3.1.3 (see also Fig. 3.4<b>). 

Fig. 3.”^ shows the variation of transitoric rate contribution 

as a function of X_ for various value of B_, B„ is fixed to 0.4. 

» TP 

For a set of Bp and B^, variation in transitoric rate contribution 
can be neglected after Xg is 10. This is due to the decreasing 
influence of (x>/T 7 in the rate equation (3.4). It can be seen 
from Fig. 3.7 » the major change in transitoric rate contribution 

can be brought by changing B^ from 0,2 to 2.0. The reason for Xg 

having negligible influence on transitoric rate contribution, 

after 10, is due to lethargic improvement in the driving force. 
Whereas increase in B^ increases the rate of impurity extraction 
and thus transitoric rate contribution increases. 

Fig. 3.8 shows the same axes but the influence of Bp . Here 

B_ is 0.6. Effect of X^, on transitoric rate contribution is same 

T S 

as that in Fig. 3.7, But increase in Bp decreases the transitoric 
rate contribution continuously. At high value of Bp, for example 
at Bp^l.S, permanent rate contribution is as high as 701., when Xg 


is 19. 





54 


Thus with one set of Bp, and Xg, it will not be wise to 
elucidate if the permanent contact reaction or transitory contact 
reaction dominates the powder injection process. Relative 
contribution is entirely depends on the magnitude of Bp, B^ and 


3. 1 . 4. £: Transitoric efficiency: 

Transitoric efficiency is the ratio of impurity concentration 
removed by the rising flux when it reaches to the top surface of 
the melt to the concentration it can remove i.c., equilibrium 
impurity concentration in the injected flux. 

Combining equation (2.8) and the transitoric rate 

contribution term in equation (2.12) we get 

B^ 


Cx 3 
= y“~ 




(3. 1> 


1 + 


with the definition of Xg (see equation 2.11) equation (3.1) can 


be rewritten as 


B. 


■» 

(x^) = Cx3i7g 


1+ ^ f 


(3.2) 


l7gCx3 in equation (3.2) is the equilibrium concentration of 
the injected flux. 


flux 



to the 


indicates how close the concentration of 


equilibrium concentration when it reaches to 


injected 


the top 
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surface of the melt. If transitoric efficiency is denoted fay E, 
then 


^ B (3,3) 

1+ ^ f 

For the pure transitoric reaction the same equation (3.3> can 
fae applied by putting f=l. 


3. 1 . 5 Opt-imum 

is significantly affected by the mass flux flox rate. In 

the section 13. 1 .21 three types of B^ values were defined. 

Optimum B^ is defined as the one just above which back transfer 

(from slag to melt) of impurity occurs. Optimum is obtained by 

putting the equilibrium concentration at t = 1, calculated from 

equation <2.29), into equation 2.22, for the known B_ and X_ 

P s 

values. Magnitude of B^ will also be affected by the Bp and Xg 
values, since shifting of concentration vs. time curve can also be 
done by these two parameters. 

Critical B_ as a function of X_ for various values of B_ is 
T S r 

represented in Fig, 3.^ , AT higher values of Xg critical B^ is 
high. In section (3.1,1) it was already discussed that the 
increase in Xg increases the difference between kinetic and 
equilibrium curve. This is the reason for the increment of 
with Xg. Because of the same reason increase in Bp decreases 
optimum B^ value. 
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Thus the identification of opticnutn value optimizes the 
powder injection refining process with respect to the variables 
like loading ratio, mass flux flow rate etc. 

3.2 Slag Carry-overs 

Usually in ladle metallurgical processes there will be 
carried-over slag from BQF on the top of the melt at the beginning 
of the injection process. The main component of BOF slag is FeO. 
For deS it is deleterious whereas for deP it is beneficial. 
Average BOF slag composition is 4e*/,CaO, 18*/, SiO^, 6V.Mgo, 4%Mno, 
7Y«Al2Q-?» 20*/iF6^Q, Prior to flux injection as there is 
some slag of different composition on the top of the molten metal, 
study of its effect, during powder injection, on refining 
processes is essential. 

When powdered flux is injected into the molten bath, on its 
rise to the top surface of the melt it will be mixed into the 
already existing top slag. As both top slag and powdered flux are 
of different composition and as there is assimilation of injected 
flux with top slag continuously through out the injection period, 
the overall top slag composition also changes. 

When equation (2,24) and equation (2.22) are compared, there 
are two additional terms i.e,, (top slag partition coefficient) 
and Y_ (Mass fraction of carried—over slag from BQF at t*0> in 

w 


equation (2.24) 



Fig. 3.10 Variation of ^ function of time 


for 


various values of carry-over slag mass. 






Soluble Al 
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3.2.1 Variation of as a function of t: 

is the overall partition coefficient of the slag which 
consists of caried-over slag and accumulated injected flux per 
unit of time. As top slag composition changes with time, r}* also 
changes. Thus the decrease or increase in affects the powder 
injection refining process. To study its variation with time two 
examples are chosen. One for sulphur and the other for 
phosphorus. 

3. 2. 1 . 1 Sulphurs 

For sulphur removal synthetic flux of composition 60% 
Ca0,-40%Al202 is chosen for injection. 17 ^ for this flux is 276 
for 0.03% soluble A1 in liquid steel. Mass of synthetic flux is 
10 Kg/<T of steel), and the mass of carried-over slag is varied 
from 10 — 50 Kg/(T of steel), the average carry-over slag 

composition is given in the section 3.2, injection rate is kept at 
1 kg/T of Steel, For the above conditions overall slag 
composition is calculated at each t (t varies from O to I ) and Cg 
is computed by optical basicity equation and equilibrium A1 
concentration of the melt is computed from the reaction 

(Fe0|+ A1 = CFe3 + <A1_0_) <3. 4) 

AS°3^4^= 4- 233.^7 j/inAe 

The calculation procedure is given in Appendix-4. In Fig. 
3.10, the variation of (v^j.) as a function of time (t) is shown. 

w O 

In this figure it is clear that the carry-over slag's partition 
coefficient is much less than that of synthetic slag. This is due 



61 


to FeQ present in the carry-over slag. Increase in < 71 *) with 

S Q 

increase in slag carry-over mass is also observed in this figure. 
This increase in can be explained with fig. 3.11 in which 

soluble A1 loss, Cg are plotted as a function of total slag mass 
i.e., mass of synthetic slag and carried-over slag. 

Computation of C_ and soluble Al in the melt is same as 
explained above, Cg increases gently along each iso Yg line, with 
carry-over slag mass, but drastic decrease in soluble Al content 
can be seen at Yg in the range 0 - lOKg/lT of Steel!, When is 

greater than 10 Kg/<T of steel) soluble Ai decreased to such low 
level that no further loss takes place. But Cg increases 

continuously. That is increase in occurs at Yg above 

tOKg/ (T of Steel ) . 

3. 2. 1 . 2 Phosphorust 

For phosphorus removal the synthetic slag composition chosen 
is SSy, CaO-23*/,FeO-22%CaF2. Mas® o'f synthetic flux, mass of 
carried-over slag and the average composition of carry-over slag 
are same as in the section 3.2. 1.1. a^ is evaluated as per the 
procedure given in appendix-4. Cp is evaluated from optical 

basicity correlation <see table 2.3>. computed from 

equation 2.55. r}p is computed to be 50. <T7p>Q ^ function of t 
is plotted in Fig, 3,lZ for various values of Yg. is 

than that of synthetic flux along and like sulphur increases 

with carry-over mass. This is just opposite to the partition 

coefficient of sulphur variation in the present of carry-over 



Fig. 3.ia 


Variation of ^ as a function of time for 

various values of carry-over slag mass. 





vs T for 
presence of 



car r y—over 




64 


slag- The increment in ('Op) is due to both C_ and a . 

” o P O 

3.3.2 Concentration vs. tliae plotss 

* 

Appearance of and in equation 2.24 forced to chose 
specific parameters, as input to modern, to discuss about slag 
carry-over effect on powder injection refining. As discussed in 
section 3.2.1 slag carry-over affects deS and deP during powder 
injection, in opposite manner. That is why, here also, they are 
separated . 

3.2. 2.1 Cl vs T plots for sulphur* 

In the fig. 3.13 <:he variation of sulphur as a function of 
time is shown for three different types of synthetic fluxes namely 
60CaQ-40Al203i 70Ca0-30Al2Q2» TOCaO-SOCaF^, in the presence and in 
the absence of carried-over slag. The selected synthetic fluxes 
have T 7 g 276, 310 and 500 respectively for 0.03% sol. A1 in liquid 
steel. In all the cases mass of synthetic flux is 10 Kg/(T of 
steel), mass of carried-over slag is 10 Kg/ <T of steel), synthetic 
flux injection rate is 1 Kg/<T of steel) and selected Bp, B^ 
values are 0.4, 0,6 respectively. Melt concentration as a 

function of time is evaluated from the equation (2.24). It can be 
seen from the figure that slag carry-over deteriorates, 
considerably, the removal of sulphur both in injection and 
stirring regions. This is due to the reversion of Sulphur from 
slag to melt. 




Degree of Des, % 


Fig. 3.15 DeS degree variation with synthetic and carry-over 


slag mass. 
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3. a. 2. 2 a vs T plots for phosphorus: 

In Fig, 3,1^ the variation of phosphorus as a function f time 
is shown for the synthetic flux of composition 
55y,CaQ-23y.Fe0-22y,CaF2 in the presence and in the absence of 
carried-Qver slag. Three values of are chosen. Calculation of 
and Cg is same as in section 3. 2. 1.2, It seems from this 
figure that carry-over slag is beneficial for deP. But if 
synthetic slags partition coefficient is greater than that of 
carried-over slag (this conditon is similar to that in deS with 
slag carry-over) rephosphorisation will occur for high value of 



Because of the reasons discussed above the negative effect of 
slag carry-over is observed, mostly, in deS by powder injection. 
3.2.3 DeS and DeP efficiency: 

Efficiency of impurity removal is defined by 

-j X 100 (3.5) 

depends as per equation (2^29) » Degree of deS is 

shown in Fig, 3.15" with total slag weighty consists of synthetic 
and carry-over slags^ as a variable. Slag system chosen is 60% 

CaO-W/fAl^O^ and initial soluble A1 of melt is 0.03%. In this 
figure. iso synthetic slag mass lines and iso carry— over slag 
masslines are shown. Along each iSo initial increase in 
with a small increase in carry— over slag mass and a drastic fall 
are observed. The initial increase is due to rise in 

sulphide capacity and the drastic decrease is due to loss of 




soluble A1 content of the melt (see Fio x it. 

rig. 0..13/. Along each iso Y 

line the continuous incre.se in degree of deS' uith incre.se in 
is obviously due to decre.se in FeO concentration in the over-.U 
slag. Similar to sulphur, degree of deP is also evaluated from 
equation (3.5) for various combinations of synthetic and 
carried-over slag. This is plotted in Fig. 3. IS. The slag system 
is 55%CaO-23y.FeO-22>lCaFj. With the exception of peat in Fig. 
3. IS, Fig. 3.16 is just opposite to Fig. 3.15, i.e. increase in 
Yj., along each iso Yg lines, increases degree of deP, but increase 
in Yg, along each iso Yg line, decreases degree of doP. This is 

due to the reason mentioned in the above section 3. 2. 2. 2. 







CHAPTER 4 
APPLICATIONS 

In this chapter the applications of the model is discussed 
with reference to the removal of sulphur^ phosphorus and nitrogen 
from the hot-metal /stee I without, and with carry-over of fixed 
amount of slag from the previous refining reactors. For this 
purpose the available information on the removal of the above 
mentioned impurities are obtained from the concentration vs time 
plots. From these concentration vs time plots the permanent and 
transitoric rate constants and their contributions towards the 
overall refining are evaluated, and discussed. About the 
organization of this chapter the rate of removal of sulphur, 
phosphorus and nitrogen are presented separately. 

4.1 DeSs 

Various research workers had investigated the removal of 
sulphur from hot-metal /stee I by injecting suitable amount of slag 
forming materials. For the purpose of the present analysis deS 
caused by injection of binary and ternary CaO based fluxes are 
chosen. The details about the type of the flux, mass flow rate of 
the flux, carrier-gas flow rate etc., are compiled from the 
available literature and given in the Table 4.1. From the 

concentration vs time plot of each investigator the values of Kp , 
f and E are determined by using equations (2.15) and (2.22). 
The calculation procedure of the above mentioned parameters is 
given in Flow-chart 2 and the source code is given in Appendix-5. 
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The values of and are derived by using equations 2.15 
and 2.22. <See flow chart-2 for details), n* for each type of 
slag is deterrnined by equation (2) of Table 2.1 and is listed in 
Table A-.l as X.g(Xg- • The value of f is calculated by 
equation <A3.5) isee appendix-3 for detailed procedure). With the 
known value of ^s* ®p are determined by using 
equation {2,14). From the known values of t. . , Kp and K_ are 
determined from the calculated values of Bp and values 
(B=kt^^^). The transitoric efficiency parameter is determined by 
equation (3.3) . The calculated values of Kp, K^, f and E are 
listed in Table 4,1. 

Fig. 4,1 and Fig. 4.2 show the variation of « vs t for the 
experimental conditions of the investigators listed in Table 4.1. 

In the literataure pressumption of Kp value or evaluation of 
Kp by dumping experiments were done to find transitoric rate 
contribution, transitoric efficiency (E) etc. But as per the 
above procedure, with the developed model, Kp value is evaluated 
directly from concentration vs time plot alone. In Table 4.1 the 
Kp values specified by respective authors are also listed. The 
comparison of evaluated Kp values with the specified Kp values 
show a match in Sr. No. 2 (see Table 4.1), but remaining two show 
a large difference. Matching of the results in Sr. No. 2 (Table 
4.1) is due to the fact that the permanent rate constant was 
specified by doing the dumping experiments under the same 
conditions as injection experiments i.e., same type of flux, same 





Model parameters evaluation from concetration vs time 
plots for Ref . CS93 
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amount, same treatment time, same gas flow rate etc. In Sr, No. 1 
the Kp value was assumed by the investigators. In Sr. No. 3 
although Kp was specified on the basis of dumping experiments but 
the conditions of powder injection and dumping were not similar. 
For the remaining experimental conditions no rate parameters were 
specified. Only the concentration vs time curves with the 
experimental conditions were available. For those also the model 
parameters i.e, Xg,K.p, K^, f and E are evaluated. 

4. Z DePs 

Similar to deS, DeP was also investigated by various research 
workers in hot-metal/steel by the injection of powdered fluxes 
under oxidizing as well as reducing conditions. The experimental 
conditions for some deP experiments are given in Table 4.2. For 
each investigators, the values of <p^ and are determined from 
their respective concentration vs. time plots and by using 
equations 2.15 and 2.22. The values of and 0 are given in 

o S 

Table 4.2. 

The variation of dimensionless concentration with time are 
shown in Fig. 4.3 and 4.4 by solid lines for CaC^ flux and 
CaQ-CaCl 2 “CaF 2 flux, are obtained by using evaluated and 0^ 
values in equation 2.22. The experimental points are superimposed 
on the respective figures. For the investigator No. 3 in Table 
4.2 the calculated values of and ©g are given, but the « vs T 
plot could not be given because of the lack of experimental points 
for comparison. For this case deP degree is calculated by the 
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following expression; 

Cx^3 - Cx^3 

deP degree = x lOO 

The calculated values of deP degree are reported in Fig, 4 5 
along with the experimental points. 

4,2,1 DeP witti CaCg flxix; 

From Fig. 4.3 it can be seen that the experimental data 
points can be reproduced well by the model equation ( 2 . 22 ) of the 
present investigaton . From their experimental result the 
coefficient of utilization of CaC^ for deP is determined. The 
coefficient is defined by 

T7gYg Y_ 

Utilization coefficient <<5) ® r — ® = r,— (4.2) 

T 7)g Y^ T 

In equation (4.2) Yj is the actual amount utilized inorder to 
decrease the concentrat ion of phosphorus from its initial to some 
final value and Y_ is actually utilized among. X_ can be 

O O 

calculated from the following mass balance equation 


m_(x,> » m. 
S f L 


Cx 3 - 
o t 


m^ 


# s 

Since (x,)* r>_Cx,3 and — « Y-, 
f S f m^ 5 


<4.3> 

after some rearrangement 


and transformation, 

„ l-cu, 

•V ^ f 

^ ‘ “ -S— 


(4.4) 


The value of ot. is determined from the concentration vs time 
f 


curve, for injection, from Ref C343. Using this valuein 

equation (4.4) X_ calculate 

S 

with phosphorus as follows 


equation ( 4 . 4 ) Xg calculated, X^ is calculated from CaC^ reactic 



Degree of Dep^ % (Measured) 



4.5 Comparison of calculated 



I L 

80 90 

Vo (Calculated) 


measured deP degree 
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SCaC^ + 2P — ♦ + 6C 

^®°4.57) * " “ 0.012T K.J/mole 


At 1879 K assuming a 


caC^ ' ^Ca P ^®°4.57) '^^lue, 


we 9^^ 

ap = 3,78xl0''^a^ 

As carbon and phosphorus concentration are low in steel, 1 
wt’/. standard state can be taken 


h„ = 3.78x 10 ^h^ 

L- =- 

From stoichometric relation of <4. 5) 
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(4.6) 


CC 3 « CC 3 + 

eq o 62 


{ CP 3 - CP 3 ^ 
o eq 


3 - LP JV (4.7> 

From the solution thermodynamics the following relations are 


valid 

h„ = f„e%P3 

r r 

he = fc«:y.c3 

r Pr 

log fp « ep C%C3 + e CXCrl 
» 0,24Cy,C3 “ 0.001Cy,Cr3 
log fe “ 0,IC%C3 - 0.036CKCr3 


(4.8) 


By solving equations (4.6), (4,7) and (4,8) simultaneously 

for the given initial carbon, phosphorus and chromium contents of 

stainless steel bath, equilibrium phosphorus (CP^^3) content in 

the melt can be determined as follows. 

In the selected example CCo3 = 0.65%, 12%, 

0.03%, and CP 3 evaluated is6«807 ppm. Using equation (4.3), in 
eq 

which « is equal to ot Ca is ohosphorousl , is evaluated to be 
T eq 
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44.5. 

In the similar way is calculated for the other initial 
phosphorus content, and these are given in Table 4.2. By equation 
(4.2) <5 is calculated and plotted in Fig. 4.6 against CaC^/T 
injected. In this figure calculated values of <5 are also shown for 
the case of dumping CaC^- It can be seen in the Fig. 4,6 that the 
utilization coefficient lay within the range 1-2.5% for dumping 
experiments, whereas this value varies from 1.8 to around 9.5% for 
injection experiments. According to kitamura et al., C611, the 

maximum <5 value for dumping is 2.7%, which coincides with the 

above calculated value for dumping. The higher & value for 

injection can be explained by considering the mechanism of deP by 
CaC^. CaC^ decomposes into Ca and C. As the B.P. (Boiling point) 
of Ca is lower than the refining temperature, its vaporisation is 
bound to occur. But the vaporisation loss is decided by the 

partial pressure of Ca at refining temperature. Injection of CaC 2 
deep into the bath decreases the partial pressure of Ca, and thus 
decreases the vaporization loss. 

One interesting result from the model is the ratio upon 

dg. It is observed that this ratio is greater than the product of 
4g and initially chosen Yg value (see Table 4.1) for 
non-vaporising fluxes and much less than the same product (see 
Table 4,2) for vaporising fluxes. If the product is designated by 
the comparison of evaluated ^g/0g and from model suggests 
whether the flux is vaporising or non-vaporising. 
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4.2.2 CaO-CaClg-CaFg flux: 

In order to evaluate the concentration of phosphorus in the 
melt as a function of time by using the model, and e values 

are required. To evaluate concentration vs time plot for the 
experimental cnditions given for deP, by using this flux, see 
T ab 1 e 4 . ^ * ^T are to be comput ed {see fl ow chart*"! ) « 

To calculate Kp and from the given experimental conditions, 
selected equations in Chapter 2 are used. To calculate X , 
precisely 7?p , regular solution model is used because of the 
inability of equations listed in Table 2.3, along with equations 


{2.55) and (2.58). v, 


p_ (see equation (2.58)) is calculated by 
2.5 


the regular solution model as follows 

RTlni^pQ = aj.^2+_pS-^ <){j,^g + \aci ^CaF^^ 

« O 2 z. 


<4. 9) 


2^- „5■^ =-250000 - 127000 llOX 
Ca -P 1 


CaF2/^CaO 




For exmaple, for 


0.8991, 0.0419, X^^p^= 0.059 


a^^ 2 +_pS+ = -358560.31 (J/mole) 
= 4.08xl0~'*^ 


The value of equilibrium constant (K) from equation (2.57) is 
100.3226. In carbon saturated iron the average value of a^ is 6 x 
10 ^ (see section 2.4, 1,1), by considering carbon and oxygen 
equilibrium in the melt. 


Thus logC- from equation (2.58) is 13.087 and T^p, from 
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equation (2.55), is 1367.4, and X_ is 32.27. 0.0236). 

The calculated value of Kp, from equation (2,77), is 0.14 

min~^. Stirring energy is evaluated from the equation (2.79). To 

calculate (3j and are needed. A^ is computed from the 

equation (2.74), fay wing equation (2.75) and (2,76) to evaluate 

t„ . The calculated value of t_ and A_ are 0-163 sec. and 3846.8 
K ri I 

2 

cm respectively. fi^ is evaluated from equation (5) of table 2.4 

along with equations (2.72) and (2.73). It is 0 -€)Sl cm/sec. 

For particle diameter of 0.1mm, usually used in powder 

injection refining, the computed is 0.0 ( 3 ^ With 

t. . is 20 min. B„and B- are 2,8 and 0.11 respectively. Similarly 
inj P T 

X_» and B_ can be evaluated for other conditions given in Table 
S r T 

4,2 for this flux, 

4.2.3 DeP hy CaO-FeO-CaFg-FegO^ fluxs 

The values of and 0^ are evaluated from the experimental 
conditions through the calculation of Bp, B^ and Xg in the same 
way as discussed above. The calculated values are given in Table 
4.2. Melt concentration at the end of injection is calculated, 
and deP degree from equation (4.1) is calculated. Calculated and 

O . ■ 

measured dsP degree are shown in Fig^ 4.5 with 45 line* T}p for 
the above system is computed from the equation (2.55) by using 
optical basicity equation (equation (2>^ Table 2*3) for Cp 
calculation. Theoretical optical basicities of oxides and 
fluorides are given in Table 2,2 and the evaluation of optical 
basicity of slag-system is explained in section 2.4,1, 
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For the given flux composition say 45-56y,CaO, 14.8«FeO, 

3 A. 62 «Fe 203 and 5 . 02 y,CaF 2 , logCp, calculated from the optical 

basicity correlation given in Table 2.2, is 7.07. toga is 

* — 

calculated to be 11.496 and y?p, calculated from equation (2.55) is 
466. With the given mass fraction of slag 0.02468, X is 
evaluated to be 11-496. Permanent rate constant is calculated 
from the equation (2.77) is 0.1222 min and transitoric rate 
constant, calculated same as the procedure given in section 4 . 2 . 2 , 
is 0.046 min”^. To calculate Bp and 8 ^, Kp and are to be 

multiplied with t (in this example it is 22 min.). But it was 

I, n j 

reported that deP proceeds then only when deSi rate is practically 
ended. In the present experimental conditions given in Table 4.2 
initial Si content of the bath is less than 0.1)4. Kawasaki et al , 
C601 reported that deSi rate depends on 

(a) rate of supply of oxygen when initial silicon is greater than 

0 , 1 %. 

(b) Silicon transportat ion rate in hot-metal when initial silicon 
is less than 0 . 1 %. 

Thus from various experimental results analysed deSi can be 
taken as first order rate reaction when initial silicon content of 
the bath is less than 0.1%. Then onwards deP is active. 
Therefore the time taken for silicon removal from initial to 0 . 01 % 
is found. It is deducted from total treatment time and thus deP 
tlegree was found. 

When initial Si content is less than 0.1%, 
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CV.Sil 

dt T 

On integration of the above equation 
Ey.Si ,3 
cy.si 3 T 


can be evaluated in the same way as giv 


en 


(4.11) 


<4. 12) 
in section 


(4.2.2). 

In the above chosen exsample initial silicon content is 0.08*/. 

and calculated K^is 0.2x34 min , and thus time taken to reduce Si 

from 0.08 to 0.01*/, from equation <4.12) is 9.31 min. Thus deP 

proceeds in the remaining time i.e,, 12,69 min. Hence B_ and B 

P T 

are 1.55 and 0.584 respectively. Final concentration of the melt 
can be calculated from the equation (2,22) as per the procedure 
shown in flow-chart 1 and deP degree is evaluated from the 
equation (4.1). Fig. 4.5 shows that the predicted deP degree is 
fairly well at i310°C and 1410°C. Slightly lower deP degree was 


predicted at 1390^C and 1375*^0 . 

4.3 DeHs 

DeN by slag/metal refining is the new concept and not 
investigated as rigorous as deS and deP. The experimental 
conditions available in the literature are restricted to only 
dumping the flux as the top of the melt. From Fig. 2.4 it can be 
seen that the magnitude of partition coefficient of deN flux is 
much lower than that of deP and deS. it was identified that 
BaO-TiO^-Al^O^ type fluxes are good for deN. The experimental 
conditions for deN, taken from ref. C383. 

To evaluate the concentration vs time for the given 
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conditions under dumping experiments, the differential equation 
( 2 . 10 ) is to be modified by putting the conditions B = o and by 
putting t = resultant differential equation is 


dCx3 

dt 


» K, 


Cx3 


Cx3 — Cx3 




By integrating the equation (4.13), we get 

l+X„ 


1 + x._ exp 


Kpt 


a 




(4.13) 


(4.14) 


As there is no equation available for calculation of nitride 

capacity of the slag system, measured value of C^a- is taken for 

the slag systems in Table 4.2, and ■ 17 ^^ is evaluated from equation 

(2.69). For values as a function of composition at ta73K 

(see Fig. 2.2 and 2.3). With the given equilibrium A1 content and 

« 

the Cj^a- values 7)^^ is evaluated. For example for the slag system 
containing lO%CaQ- 35 */.BaO- 45 XAl 2 Q 2 ~ 10 >lTi 02 , is calculated to be 
72.05 with the equilibrium A1 content of the melt 0.15%. In 
induction stirring /3p varies from 0.02-0.03 Cm/sec C3B3. 
Slag/metal interfacial area, in induction stirring, is the reactor 
cross-sectional area. Taking is 0.02 Cm/sec and interfacial 
area, for the above example is calculated from the crucible inner 

*Tt "T 

dia, 06.6 Cm'^ and V is 1282,03 Cm (for 10 Kg melt), Kp is 

fi ^ < 

calculated to be 0.081 mm . Dimensionless concentration of the 
melt at time t is evaluated from equation (4.14) and shown as 
solid line in Fig. 4.7 with the experimental points. The close 
approximation of calculated and experimental points suggest that 
the model prediction is fairly well. 
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CHAPTER 5 

CONCLUSIONS m) SUGGESTIONS FOR FUTURE WORK 

AS a result of the present work, the following conclusions 
are made: 

j Matheniat 1 ca.i model , developed in the present work, can be 
used as a tool to study the influence of powder infection 
parameters in impurity extraction from the melt, 
g, Transitoric rate contribution in the overall refining process 
depends on the magnitude of Bp, and 

3 , The model demonstrates the deleterious effect of slag 
carry-over on powder injection refining process, especially 
with respect to deS. 

4. Model is the mediator between powder injection refining 
parameters and melt concentration as a function of time, 
i.e. , by using concentration vs. time curves powder injection 
refining parameters can be obtained or vice versa. 

From the analysis of the model and the equations available in 
literature to calculate model parameters the following suggestions 
for further work are made: 

1- Thorough analysis of critical stirring energy required for 
the emulsification of slag into melt, emulsified mass, drop 
diameter as a function of stirring energy are needed to 

evaluate Kp. 

Investigation on injected powder - melt contact area is 
desirable. 


2. 
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Table 3.1 r Sulphide Capacity Correlations 

corre 1 at ion | Remarks 



C39! 105 C^. 3.'t4<Nj,^j,«.IN„^g-0.8N^ 

-''s.o ^ 


C403 loc Cg=-13,913+42.B4A-23.82A^- 
-O . 02223310 -O . 02275 A 1_Q^ 

CA<0.8) 

C^= -0.6261+0.4e08A+0.7197A^ 

+0,0005144F6Q 

<A>0,8) 

C413 log Co= ^2,82-13300/T 

wh€r-9= ” 

B=5 .62<*/,CaQ)+4. 15 (y.MgO)-t . 15<y.SiQ2) 
+1 .46<%Al203> 


2.59 Fe-0. 12-0. 15% 

S alloy was 
equilibrated with 
7 q slag of calcium 
aluminate base in 
a magnesite 
crucible at 1575, 

1600 and 1650^^0, 

Si 0^=20%; 

Ca0=38~55%, 

h. 45% and 

Mg0=5%, 

The equations are derived 
from large number of 
experimental data of wide 
range of CaO based slag 
systems by using the optical 
basicity concept. 


The equation was derived by 
using the correlation 
established by Duffy et al 
C421 and the data of 
Venkatadri and Bell C433, 
Abraham and Richardson £443, 
and Turkdogan £453. 


A= (%CaO) + l ,39 (%H 90 ) + 1 ,87<ytSi02) 

H-1 ,65 <%A1 0 ) 

log Cg=lag - 1,375 






CaF 




Table 2*3: Phosphorxis Capacity Correlations 


C463 RTlni^^ 



correiat ion 


-7500X. 


+3500X, 


•4L X « W 

60000X^ --24500X- -X^ ^ + 

FeO CaO FeO MgQ 

22500X_ n -33600X^ _ n “ 

FeO SiO^ FeO^ ^ MgO 

4800X^ ^ X^ ^+15700X- ^ ^ - 

Fe0^^3 (I 9 O FeQ^^3 SxO, 


lag Cp=-18, 184-+35,84A-22,35A + 


229 30 A 
T 


- 0.06257Fe0“0.04256Mn0+ 
^0,3 


O ,3598^0^ 


P {%P) 

P Q 


I o, Cp .0 . 5 1 < 23N(. 1 ™„ga*eNp . 

2 2 5 

2V920/T-19.2a+la9C(XP)/Np^^ 3 

log = ,0.478 lag N + 

Np o 

2^5 

1«955 far CaO-FeQ-SiQ^-P^O^- 
<CaQ/SiQ^=3 F6^Q=20y.) 

=^0.4791ag ^ +1.813 for CaQ- 

Fe^Q-CaF2-P205 

(CaO/CaF«=:l; Fe.Q=20%) 

r ^ tl-i 


c =— 




Remarks 


P+2.5g=(PQ2 


log Kp =17060/1-8. 51 
Slags System : 

FeO-Fe^O^-P^O^-SiO^^CaO-Mg^ 
<from 1380-16a0°C) 


Derived from literature 
data and then tested by 
performing experiments. 
Recommended optical 
basicity values are to be 
used . 


7-13%Na2Q & 4XBaO 
additions are done. 

Experiments are done at 
1550°C 

Using AG for the followin 
equation 2P+5g = (P^Q ^) 5 

the following 

relationship is derived, 

- _ AG° 1 . , 

*^9=- 2RT - 2^^°^ 

- ^ ^ w 


and lag v 


is taken from Turkdogan's 
treatment • 





Tikblo 4x Mjiss t,r iirisf er coof fie. lent, reltttloris for drops 







Table 2. 5t Slag-Systems and the evaluated rate constants 


Sr. 

No. 

Slag System 

Kpimin S 

15 

657 . 86203 - 2 , 5y.CaO- 10%A 1 2 Q 3 
(Inductively stirred) 



— *■ Sulphur 

0 . 0273 


— > Phosphorus 

0 . 06634 

16 

757 , 88203 - 257.030 
(Inductively stirred) 



— > Phosphorus 

0 . 033 

17 

m 

Na2C02 <1400°C) <£=76W/T) 



— * Sulphur CSi^=0.0173 

0,2283 


— * Sulphur CSi^=0.133 

0.2283 


— » Phosphorus C Si ^=0.0 17. 3 

0.0691 


— *■ Phosphorus CSi^=0. 137.3 

0.05066 

18 

m. 

Na 2 CQ 2 (1300^0 <^^=76W/T) 



— ♦ Sulphur CSi =0.017.3 

Q 

0.0723 


— >• Phosphorus CSi^=0.0173 

0.0146 

19 

1 

3 1 7,C aQ-337.S i Q 2 - 1 77, A 1 2 Q 3 - 1 97.0 aC 1 2 

1 

(Al-deoxidized 1350 C) 

0.2533 

20 

Na 2 C 0 ^ CPellets3 

0.2533 


CPowde r 3 

0.194 

21 

(40.607,Ca0)- ( 15 - 307 .Ca 82 >- 

0,035 


(25-407. (8e0) 

— > (deP) 

(250T arc process) 

0 . 084 






Table 2.5s Slag -Systems and the evaluated rate constants 


Sr . 
Nc. 

Slag System 

Kp <min 

22* 

6T,V- =1.2, B=1.73 

Ar 

0.39 


6T,V. =1.8, B=l,96 
^Ar 

0.55 


5T.V, =0.9, B=2.43 

Ar 

0.45 


* 

5T.V^ =1.1, B=2.75 

Ar 

0.44 


‘ 4.5T.V^ =^0.075, B = ? 

Ar 

0.02B 


(Si-^deox> 



4,5T.V^ =0.075, B=? 

Ar 

0.042 


(Al-deox) 



ISOT.V^ = ?, B= 3 

Ar 

0.066 


<Al-deQx) 


23 

aoy.N a^C0_- 1 3 , 37.3 i Q--6 . 67y,Fe„Q_ 

2 3 2 2 3 

0.5 





values 

. 3 
dm n 


y 


are 


taken from the ref. 


' 



'^CaO * ^CaF 




V 


in 


t . sec 


B 








Table a. 6s Evaluated parameters for hot metal experimental corulitions 


Sr, 

No- 

s 

LoadiriQ 

Ratio 

Y 



Stafee of 
at hot metal 






temperature 


(W/T) 

3 

(KQ/tT\ > 



<min ^ ) 

<1623K-1673K) 

1 i 

b.22 

666.67 

0.0067 

1060 

0.245 

Sol id 




0.01-0.02 

542“ 



2 




1-35 



3 : 

248.77 

15-55 

0.0058 

CS 3 = 
eq 

,05ppm 

.0374 “ 

, 05 6 


4 

35 

78.8- 

. 006 1 4 - 

405,5 

.068 “ 



121.2 

.0095 


. 1044 


5 

255 

10-40 

.02-. 05 

25 

,0257- 
, 1026 

Liquid 

6 

12.25- 

128.57- 

,0029“ 

1507 

.0734- 



17.14 

300 

.0051 

1 

1 

1167 

. 1458 
,0734- 
. 1458 


1 

1 

: 1 

f 1 

1 

1336,6 

.0734- 
, 1458 

Solid 

7 ' 

120“ ! 

3 , 33“ 

.0018 


.006 - 



180 

12,5 



.0183 


8 

400 

19,725 

.005 

CS 1 = 
eq 

,0612 






,05ppfn 

1 





21 .375 

,041 


,0662 

' Liquid 

9 

950 

4,17 

,0294 

CS 3 = 
eq 

,0155 






.05ppr;^ 




340 

10,42 

,0294 

CS 3 = 
eq 

,025 






.05ppm 




1130 




,0104“ 

,02 



1130 



177 

,0178“ 

,03 

L iquid 

10 



.0332 

177 


L iquid 




.0051 

CS 3 = 
eq 

, 05ppm 


L iquid 


Conttd 


Table 2-65 Evaluated parameters for hot metal experimental conditions 


Sr . 
No • 

£ 

(W/T) 

' "'*1 

Loading 

Ratio 

! <Kg/m^) 

i 

Y 

* 

V? 1 

1 

’ 1 

(min 

State of slag 
at hot metal 
temperature 

<1623K“1673K) 

11 



.02 

CS 3 = 
eq 

, 05ppm 

i 


12 

41.4 

100 


,862 

^ .0876 

Solid 

16 


450 


tS 1 = 
eq 

.07ppm 

f 

i • 


17 

1165,6- 

1795,66 

2-3,2 



^ .01244- 

i .01612 

1 

Solid 


600 

6 

.00113 


1 

,02183 

Solid 

18 

976,35- 
1285 ' 
274 

10,52 

5.76 

.0428- 

.1045 

.06 i 

500- 

1100 

.02542- 

.0717 

.0252 

Liquid 

Liquid 


^4ate: Sr. No. in this Table is same as that in Table 1.1 











Table 2.7: Evaluated p^rsun&terrs for Steel experimental coridltlonsi 


Sr. 

No. 


Loading 

Ratio 

V 

* 

n 


Stats af slag 
at stssi making 







temp 3 rat are 


(W/T> 

CKg/m^) 




(1873^K) 

1 

121.42 

39,13 

,0048 

1555,47 

,06474 

Partially Liquid 

2 

36,77 

65,75- 

,0028 


,061- 

Partially Liquid 



78,125 



,0724 




69,375 

,0028 


, 06428 

Partially Liquid 



62,5- 

,0028 


.058- 

Partially Liquid 



64,375 



.06 




57,375 

,0028 


,05316 

Partially Liquid 



55.875- 

,002- 

CS 3 = 

,052- 




71 

,003 

eq 

,23ppm 

,066 


3 



,015 

113,82 


Liquid 





2103,35 


Sol id 


' 1 



131.2 


Liquid 

4 

53.6- ! 

73,86- I 

,013- 

575 

,0975- 

Liquid + CaO 


78,6 ' 

125 ; 

,015 


, 137 


5 

! 

264 

21 

,0042 

1 300 

,0384- 

.0577 

Sol Id 

6 : 


,77- 

,02- 

300 

,01 1- 

Liquid 



1,54 

,04 


,022 


7 

106 

60— 

,03- 

46.13 

, 0082- 

Liquid <40-50% CaO) 



200 

, 1 


,0272 

Liquid + CaO 050%) 




,001- ' 

46,13 


Liquid (40-50% CaO) 




,005 



Liquid + CaO 050%) 

8 



,01- 

,04 

CP 3 = 
eq 

8,2ppm 



10 

42-63 

37-42 

,0015- 

CS 3 = 

,04- 





,0018 

eq 

,23ppm 

,052 j 






Table 2. 7i Evaluated paraineters tar Steel experimental conditions 


Sr, 

- 

Loading 

Y 



State of slag 

^4o• 

£ 

> Ratio 


i ^ 

1 

at steel making 
temoe nature 



1 <W/T) 

(Kg/m^) 



(min ^ ) 

{1S73°K) 

11 

i 


1 .015- 

■CP 3 = 


1 



\ 

,045 

eq 

8ppm 

i 

1 

12 

29,35- 

20-36 

,0008- 

CS 3 = 

i 

.0313- 

1 

i 

1 

1 

88,94 


,0011 

eq 

^ ,23ppm 

,03 

t 

1 

1 

32,43 

36-40 

,00062- 


, 14- 

i 

l 




,00154 


,35 

f 

t 


43,4- 

20-23 

.0004- 


,021- 

Lieu id 


115,42 


,0026 ! 


.0374 



Note: Sr, No. in this labia is same as that in Table 1,2 













CtlcuUted 











Table 4.^1 DeP Evaluation 














APPENDIX - t 

General solution for powder injection refining without slag 
carry-overs 

Governing differential equation for concentation change of 
impurity in the melt with time during the powder injection when 
slag is not carried— over is (see equation (2.15)) 


d(«-l ) . , 1 V fa . . 

-3i — I - ^*s 

It was already discussed in the 
of equation <2. 15) is similar to 
equation . 


<2.15> 

section (2.2.1) that the form 
that of Leibnitz's linear 


dY 

dx 


+ Yf<x) = C 


( 2 . 21 ) 


The general solution of equation (2.21) is given by 


Y = e 


-Jf(x)dx;c dx + C^l 


(Al.l) 


where is integration constant. 

By putting similar term from equation (2.15) into equation 
(Al.l), we get 


-f) dT f 

(«-l) = e < (- 






f 








9s 


dT+C 


.} 


dT 


=1' 


-{^gT T^glnr) 


-( 4 fr_T +0 Inr) f ^„T e. 




s s . 

e T dr 


1 + 2 




dr 


(0g+l) 

00 . . s ® ^ < ^e'*’ ^ ^ 

<^„t) T S 


2! (Bg+S) 


- 2L < 


m=0 


l+0g+m) 


(A1.2) 


(A1.3) 


(A1.4) 




i , e . , = O 

Thus, 

“ /J. -r \*** 

Z <^gT> 

<n'. <iH-©:'-fii — 

m=0 “ 

Equation (2.22) is the obtained general solution. 


(A1.5) 


< 2 . 22 ) 



APPENDIX ” 2 


Numerical solution for powder injection refining withwn^ slag 


carry-over . 


Boverning differential equations for concentration change of 
impurity with time during powder injection refining with slag 


carry-over from BOF is (see equation 3.20) 


d (ct- 
dT 


“ = {v ^ ^ ]} 

/?s L ® J 


<«— i ) 


(2.20) 


By using crank-Nicholson 's FDM 

'"j+r "j’ _ f. . ®s 

■ P Piilli* 




u 

fr . . . ■‘•T 
3+i 3 

2 


f "I 

j + 1 J . 

2 


(A2. 1 ) 


Defining (r > = (t..,+ t.) 

^ av j+1 j+1 j 


ir*) 

Q a 


= (7?*) + irt*) . 

V '03 + 1 0 3 

Ji+1 L i 


(2.25) 


ot .-a . , . 

3 3 + 1 


. Sg ^ 


- > 1 - 


After rearrangement 

-l)Cl+T3 = <«.-l)Ci-T3 - <^ct^T 
3 + 1 3 S 

T is defined in equation (2.15) 

Equation (2.24) is the obtained numerical solution. 


(A2.2> 


(2.24) 



appendix 3 


Transitoric Rate Contribution CO 

f is rate contribution of injected flux in the overa 

refinin. process. Therefore f is the ratio of area under 

transitoric rate curve to area under total rate curve. 


(Rate^) dr 


T = l 

f<Rate )dT 


(A3.1 ) 


Considering 
reqritten as 


-j ^ i mo interval dT. ^ A 3 • 1 ) c an be 

idering the equals time interval ui , 


<Ra^®T^T+dr 


CA3,2> 


(Rate > 


T+dT 


From equation <2,15> 


Rate = <oi-l) <■^>5 + “ 


From equation (2.12) 
Rate^ = ~ 

Then , 


(A3. 3) 


(A3. 4) 


Vo (.^Q + 

f feb (A3. 5) 

« der^ds as ^3, <^3 and ^3 depend on B, , B,, X3 and f . 

Therefore f depends on only Bp, B^ and Xg. 

(A3. 6) 

Thus f = fliBp, B^, Xg) 


(A3. 6) 



APPENDIX 4 


Calculation procedure of a^: 


If the slag in contact with the melt contains 
the components CaO, MgO, FeO, MnO, SiOg, AlgOg, then the mole 

fraction of oxygen ions of the slag melt is determined as follows, 


" "CaO^'^MgO'^^Fea'^^hgO ^"siOg ^AlgOg ^"PgOg 
n 


(A4. 1 ) 


"SiQg ♦ ^"Alg 03 " 


n 


(A4.a) 


2- 


Relation between X g- and activity coefficient of FeO is C593 

o 




S- 


1A75 

T~ 


( AA , 3 ) 


<b 1 tj U 


1 og 


iFeO) 


+ E.60A 


( AA . A ) 


Thus activity of oxygen in melt can be determined 
equation (AA.A) for the given slag composition of the slag. 

E . T . Tur kdogan 111] suggested the following solubility p 
relation for the Al-0 equilibrium for AJ^ killed steels 

CAll'^CO:^ = l.Sxl0“^^ 

Final content of the melt when the initial A_l_ ( A1 

initial Q ^ 5 ace given, is 

5A 

AJL. = 81. - 

f o AS 


□ 

•no 


- A.93SAX10 


-5 


f rome 


roduct 


( AA . 5 ) 
) and 


(AA.6) 



while (factor >= 0.0001) 

div = factl — 0.0; 
if (m > 0) 

for (k = 1.0; k <= m; k++) 

factl += log(k)-5.0; 

(j i V +“ 5 - 0 ; 

SLim *+*” "f actor 5 

fn+‘+* ; 

fact = 1 - 0 ; 

xCilCpointl = l. 0 -phi*t*exp(-phi*t)*sum; 

fra = sum_rate_tra/sutn_rate; 
printfC'fra = */.f \n" , f ra ) ; 

number = point; 

^ printf'phi . V.fNtthPt3 . «fXn“.phi,th;taM^ 

for (t=1.01,paint=number;t<-max_time, 

du».yi = 

dummyS - ' ■ (-dimmyl ) ) / ( 1 . 0+laP’daC 1 3 > i 

X C i 3 Cpoint 3 , r ^ i*-i- i n i* ( xC i 3 ^ 

rate_perC i 3 Cpoint - P ^ ^ 0 +lamdaL i 3 >-l • 0 > ^ 

(1.0+lamdai:i3)-1.0)/lamdaCi3, 


rateCi3Cpoint3 = rate_per!:i3i:point3 ; 

} 

> 

if (equ_super == ’ y ‘ > 

i 

no_lamdas =1; 

IvCno curves! - lamdaLv , 

if ^nd .curved >1) (na.curves - 2)1 

i 

i = 0; V 1 =^rt=.ri + 13 ) 11 (lamdaCi3 < 

no_lamdas = ^ ^ i + 1 3 > ) ? <no_lamda 5 +l ) -.no.lamdas; 

.1 _ f/iamriaCi3 > lamdaCi+13) 11 

lytno_curve5+no_lamda5^13^^^^^^^._^^^^^ ^ lamdaC i + 1 3 : lamdaC 

else if (no.curves >2) 

(i . Oi i <= no-puf'-es 


- S; i++) 



i 

no lamdas 


(damdaCi] > lamdaCi + 13) II (lamdaCi: < 

lamda[i+13 ) ) ? (no_lamdas+l) :no_lamdas; 


IvCno curves+no lamdas-d = ((lamdaCiH 
^ ' < lamdaCi + n) ) ? 


> lamdaCi + n) I I (lamdaCi^ 
lamdaC i+1 3 i lamdaC i 3 ; 


3 

> 

else if (equ_super == 'n') 
no lamdas = 0; 


ifp = fopen("KINETIC_DATA","W); 

no.points = fabs(max_time/time_interval>+l; 

if ( type_of_graph == 1) 

fprintf ( ifp , "•/♦i\n" ,no_curves+no_lamdas) ; 

for (i = 0; i <= na_curves - 1; 


i 

fprintf ( ifp , "%i\n“ , 
fprintf (ifp, "•/♦i\n“, 

fprintf ( ifp , "LY=y.3. 


i + 1 ) ; 

na_points ) ; 

af_Bp=y.5.^f_Bt=%5. 4f\n*' 


lamdaCiljkpCil-M-tinj , 

ktCi3*tinj ) ; 


for (t = 0.0, j = 0; t <= max_time; t +- tifne_interval, j++) 
fprintf ( ifp , "3* f \ t Xf \n " ,t,xCi3Cj3) j 
> 

for (i = na_curves; i <= no_curves+no_lamdas - 1; i++) 

{ 

fprintf ( ifp, "*/.i\n", i + 1) ; 
fprintf ( ifp , "%i\r\'' ,no_points) ; 

fprintf ( ifp, "LY=7,3.Ef\n",lyt: i3) ; 

for (t = 0.0; t <= max_time; t +- time_interval) 
if C t < 1 .OOB) . - ^ V , 

fprintf ( ifp, "%f\t'/.f\n", t, 1.0/ (l-0+lyCi3*t) ) ; 

else i"F > 1-008) 

fprintfC ifp,"%f\t%f\n" ,t,1.0/(1.0+lyCi3) ) ; 

} 

griph^set t ings ( max_t ime , no_curves , no.lamdas , type_ot _graph ) ; 

> 

else if ( type_of_graph == 2) 

fprintf (ifp, "yi\n", 3 *nD_curves ) ; 
for (i =0; i <= no_curves - Ij 
{ 

fprintf (ifp, "’/.i\n", i + D; 

fprintf ( ifp . .RATE.TDTADin" , 

fprintfllfp, LY-/.3.2f_Bp [ i j.t in j , kt C U.t inj ) ; 

for (t=time.intpryal,i=l;t<=i"a« ti.nEjtt=time_interval, ir+> 

fprintf ( if p , '"Af\VAf\u" , t ,rateC 1 3 C j 3 ) ; 

> 



1; i++) 


for <i = 0; i <= na_curves — 

<. 

fprintf<ifp, "%i\n“, no_curves+i+l ) ; 
fprintfCifp, "Xi\n", na_points-l ) ; 

f pr i ntf ( i f p , " L Y=X3 , 8f _Bp=*/5 . _B t=y,5 . _ ( RATE_P ) \n ” , 1 amda C i 1 , 

l<pCin*tinj,kt[i]*tinj) ; 

for ( t=time_ interval , j = l ; t<=max_t ime ; t+=tinne_ interval , j++ ) 
fpr intf ( if p , "%f \tJJf \n" , t , rate_per C i 3 C j 1 ) ; 

> 

no_paints = f abs ( 1 . 0/t ime_interval ) ; 
for (i = 0; i <= no_curves — 1; i++) 
i 

fprintf(ifp, "*/.i\n", 2*na_curves+i + l ) ; 
fprintf(ifp, "*/li\n‘', na_points); 

fpr intf ( i f p , " LY=*/.3 . Sf _Bp=%5 . 4f _B t=7,5 . iff _ ( RATE_T ) \n " , 1 amda C i 1 , 

kpCi3*tinj,ktCi!l*tinj) ; 


for ( t=t ime_interval , j=l j t<=l . 0 ; t+=t ime_interval , j++) 

fpr intf ( if p , "7.f \t7.f \n'' , t ,rate_tra[I i 3 C j 3 ) ; 

> 

fclose( ifp ) y 

graph_sett ings ( max_t ime , no_curves , no_lamdas , type_of_graph ) ; 

> 

5ystem( "gr " ) ; 

> 

i nput ( ) 

C 

int i ; 


system < " c 1 ear " ) ; 
printf ( "\33h\33J" ) ; 


printf ("\33&,aSr6C DIMENSIONLESS TIME RANGE = 0.0 - 
scanf ( "*/.f " , &!max_t ime ) ; 

printf < "\33&:a4r5C "); 

printf ( "\33&a7r5C 

printf ( "\33&a5r6C 1. CX3/CXo3 Vs t* "); 
printf ( "\33Sca6r6C 2. RATE Vs t* " ) ; 

printf ( "\33S<a3r6C TYPE OF GRAPH ? "); 
scanf ( "7.i " , &ctype_of _graph ) ; 

if ( type_of _graph == 1) 

C 

printf ( ■■\33&a8r6CNUMBER OF CURVES ( EXCLUDING EQUILIBRIUM CURVES)?" 
scanf ( "%i" ,&tno_curves ) ; 

> 

else if ( type_af _graph == 2) 

{ 

printf ("\33&a8r6C NUMBER OF TOTAL RATE CURVES ? "); 
scanf ( "%i" ,&no_curves) ; 

> 

printf ( "\n\n\n" ) ; 



for (i ■ Oj i <■ na_curve» - 1| i++) 

i 

printf ( "\tCURVE NUMBER ; y,i\n" , i + 1 ) ; 
printf ("\t\tLAMDACya] =",i+l); 
scanf ( "Vif " j&lamdaC i ] ) ; 

printf("\t\tKpC%i: ( /min ) = 
scanf ( "'/if " , S>kp C i 3 ) ; 

pr intf ( " \t\tKt C'/i i 3 ( /min ) = + 

scant ( "'/»f " ,&<kti:i3 ) ; 
printf < "\n" ) ; 

> 

printf ("\n\tINJECTION TIME (min) ? "); 
scant ( "'/.f " j&tinj ) ; 

getchar ( ) ; 

if ( type_af_graph == 1) 

< 

printf ("DO YOU WANT TO SUPERIMPOSE THE EQUILIBRIUM CURVES (y/n) ? 
scant ( "VtC " , &cequ_super ) ; 

> 

> 

graph_set t ings ( max_t i me , no_curves , nD_lamdas , type_af _graph ) 
float max_time; 

i nt no_curve5 , no_ 1 amdas , type_af _graph ; 

#define na_ticx_inter 10 
ttdefine no_ticy_inter 10 
C 

FILE *df p , *if p , *f fp , *gf p , *fopen ( ) , *1 tp , *lfp , *mfp , *txtc , *t i tf , »out 
int i , 1 type , m , no_cr ; 

dtp = fopen( "DATAFILE" , "w" ) ; 
fpr intf (dtp, "KINETIC_DATA\n" ) ; 
f close ( dtp ) ; 

ifp = fopen( "XYRANGE" , "w" ) ; 
ftp = fopen( "TIC_INTERVAL" ,"w" ) ; 
if ( type_of _graph == 1 ) 

{ 

f printf ( if p , "’/*f \t 1 .0\n0. 0\t0 . 0\nm" ,max_t ime ) ; 

fprintf (ftp , "'/if \n%f \n" ,max_time/no_ticx_inter , 1 .0/na_ticy_inter ) 

> 

else if ( type_of _graph == S) 

C 

fprintf ( ifp , "'/.f\te.0\n0.0\t0.0\nm" ,max_time) ; 

fprintf ( ftp , "'/.f \n'/.f \n" , max_t ime/no_t icx_inter , S . 0/no_t icy_ inter > 

} 

fclose ( ifp ) ; 
f close (ftp ) ; 

gfp = fopen("XYLABELS","w ")5 
f pr intf ( gfp , "*/**/*£ . If \n%%2 .lf\n"); 
fclosB(gfp); 



1-Fp = fopen("LINECOLQR","w") ; 
mfp = fopen ( "LINEMARKER" , "w" ) ; 

Itp = fopen ("LINETYPE","w"); 

m = E;ltype = 0; 

fprintf ( ifp , "*/.i\ni' ,m) ; 

if ( type_af_graph == 1) 

no_cr = no_cur ve5+no_lamdas ; 

else if ( type_of_graph == 2) 

no_cr = 3*na_curve5; 

for (i =0; i <= no_cr-l; i++) 

{ 

fprintfdtp, "y,i\n" , Itype) ; 
fprintf (mfp , "0\n" ) ; 
if ( m < 6 ) 

else if (<m == 6) II (m > 6)) 
m = 2 ; 

if (Itype < 6) 

Itype += 2; 

else if (Itype > 5) 

Itype = (Itype == 6) ? 1:0; 
fprintf ( ifp , "%i\n'‘ ,m) ; 

> 

f close ( If p ) ; 
fc lose ( mfp ) ; 
fclose( Itp ) ; 

titf = fapen( ’’TITLES", "w" ) ; 
if ( type_Qf _graph == 1) 
fprintf ( titf , ■’t*\nCX:/CXa3\n" ) ; 
else if ( type_of _graph == 2) 
fprintf (titf , '’t*\nRATE\n’' ) ; 
f c lose (titf) ; 

txtc = fopen ( ’’TEXTCOLOR’’ ; 
f pr intf( txtc , ’’ 4\nA-\n" ) ; 
f c lose < txtc ) ; 

outf = fopen ( ’'0UTPUT_CH0ICE’' , "w" ) 
fprintf (outf, ’’screenXn’’ ) ; 
fclose(autf ) ; 
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